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ABSTRACT 

We use stellar population synthesis modeling to analyze the host galaxy properties of a sample of 
33 UV-selected, narrow-lined active galactic nuclei (AGNs) at z ~ 2 — 3. In order to quantify the 
contribution of AGN emission to host galaxy broadband spectral energy distributions (SEDs), we use 
the subsample of 11 AGNs with photometric coverage spanning from rest-frame UV through near-IR 
wavelengths. Modeling the SEDs of these objects with a linear combination of stellar population and 
AGN templates, we infer the effect of the AGN on derived stellar population parameters. We also 
estimate the typical bias in derived stellar populations for AGNs lacking rest-frame near-IR wavelength 
coverage, and develop a method for inferring the true host galaxy properties. We compare AGN host 
galaxy properties to those of a sample of UV-selected, star-forming non-AGNs in the same redshift 
range, including a subsample carefully matched in stellar mass. Although the AGNs have higher 
masses and SFRs than the full non-active sample, their stellar population properties are consistent 
with those of the mass-selected sample, suggesting that the presence of an AGN is not connected 
with the cessation of star- formation activity in star- forming galaxies at z ~ 2 — 3. We suggest that 
a correlation between Mbh and galaxy stellar mass is already in place at this epoch. Assuming a 
roughly constant Eddington ratio for AGNs at all stellar masses, we are unable to detect the AGNs 
in low-mass galaxies because they are simply too faint. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: 
active galactic nuclei 



1. INTRODUCTION 

Active galactic nuclei (AGNs) may play an important 
role in the evolution of their host galaxies, as evidenced 
by the the tight correlation betwe en black hole mass 
and bulge mass in local galaxies (iFerra rcsc & Merrittl 
[2000tlGebhardt et al.ll200lIIGultekin et al.ll2009D . AGNs 
and star formation are both fueled by gas accretion pro- 
cesses, and the gas content and temperature in galaxies 
can be influenced by the energy output from an AGN. 
This latter process is often used as an explanation for 
the distribution of luminosities and colors among mas- 

1 Based, in part, on data obtained at the W.M. Keck Obser- 
vatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California, 
and NASA, and was made possible by the generous financial sup- 
port of the W.M. Keck Foundation. 
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sive galaxies in the local universe. Feedback from an 
AGN prevents further star formation, causing galaxy 
colors to redden as well as contributing to the expo- 
nential fall-off at the bright end of the galaxy lumi- 
nosity function (iHopkins et alJl2006b iCroton et alJl2006t 
iKhalatvan et al.l 120081: iSomerville et al.l |200§T How- 
ever, strong evidence linking the presence of an AGN 
to the cessation of star formation is still missing (but 
see lTremonti et af1l2007t ). One popular method that has 
been used to examine this possible connection is the com- 
parison of AGN host galaxy stellar populations to those 
of similar, non-active galaxies. With a controlled ex- 
amination of the stellar populations of galaxies hosting 
AGNs, we can understand both the influence of an AGN 
on galaxy properties as well as the timescales involved in 
AGN activity. 
AGN host galaxy properties have been studied across 
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a wide range of redshifts. In the local universe, 
results from AGNs identified based on their opti- 
cal properties indicate that these objects are hosted 
by massive (\og(M*/M p,) > 10), early-type galaxies 
(jKauffmann et all l2"b03) . Similar trends are observed for 
local radio-loud AGNs; the radio-loud AGN fraction is 
obser ved to be > 30% for galaxies with log (M*/M Q ) > 
11.7 (|Best et all 120051 ). Studies of optically- and X- 
ray-selected AGNs at intermediate redshifts (z ~ 0.5 — 
1.5) indicate that these AGN host galaxies are mas- 
sive, with characteristic ste llar masses of logfM^/M^ 
~ 11 (TBundv et al.l [200l lAlonso-Herrero et all I200S 
iBrusa et al.l 12009ft . If AGNs are responsible for quench- 
ing star formation, this effect should be reflected in the 
optical colors of the host galaxies. Optically-selected 
AGN hosts in the local universe are foun d to lie on 
the red sequence (jSchawinski et all l2007| ). while, at 
higher redshifts, both optically-selected and X-ray se- 
lected AGNs reside in galaxies in the transition re- 
gion between the blue a n d red sequences dNandra et al 
20071; iSalim et al.l 120071: iCoil et al l 120091; Ifflckox et al 



2009t JSchawinski et al.ll2010h . However. iSilverman et al 



(2009) found that X-ray selected AGNs at z ~ 1 often 
reside in galaxies with higher star-formation rates than 
similar-mass non-active galaxies, and that the AGN frac- 
tion is higher in the blue cloud. These results highlight 
the importance of comparing AGN hosts with control 
non-AGN samples matched in stellar mass in order to 
mitigate the selection effects associated with missing low- 
mass black holes in small galaxies. 

To uncover evidence of ongoing processes linking AGNs 
to their hosts, it is important to examine AGNs at even 
higher redshifts, during the epoch when the bulge stellar 
population was forming. At z ~ 2 — 3, both the star- 
formation-rate density and the black hole accretion rate 
peak in the universe, and this era is ideal for examining 
the relationship betw een AGN host gal a xies and simi- 
lar n on-active hosts (|Reddv et al.l [2008; Hopkin s~et al.l 
2007). At these redshifts, it has proven very challenging 
to derive robust constr aints on the host-galaxy properties 
of un obscured AGNs (|Targett et al.ll2012t ISantini et al.1 
120 121 ) . Therefore it is preferable to limit the study 
of host-galaxy properties to obscured AGNs, in which 
the emission from the central engine does not com- 
pletely outshine t he light from the stellar population 
(|Assef et al.l 12010ft. Results at z ~ 2 - 3 for UV- 
selecte d AGNs (lErb et all l2006bt ) and near-IR selected 
AGNs (jKriek et al.ll20W indicate that narrow-emission- 
line selected AGNs at these redshifts are mainly found 
in high mass galaxies (log(M*/M Q ) ~ 11). These re- 
sults are similar to those for X-ray selected AGNs from 
both the Cosmic Assembly Near-in frared Deep Extra- 
galactic Legacy Survey (CANDELS; iGrogin et aTll2011b 
IKoekemoer et^muh . and the GOODS and COSMOS 
fields, which also provide evidence that A GNs are hosted 
by massive (log(M*/M^) > 10) galaxies (iMainieri et al.l 
[20111 Kocevski et "all 12011 ISantini et al.l 12012ft . These 
authors also estimate host galaxy colors for AGNs that 
put t hem both i n the blue cloud and the transition re- 
gion. iXue et ahl (|2010l ) carefully examined X-ray selected 
AGNs out to z ~ 2 — 3 and found that, when compared to 
a mass- matched sample of non-active galaxies, AGNs and 
non-AGNs occupied similar regions of color space, be- 
coming bluer as redshift increases. These results demon- 



strate that host galaxy properties for large samples of 
AGNs need to be carefully compared to the properties of 
a controlled sample of non-active, star-forming galaxies 
to understand the effects of an AGN on galaxy evolution. 

In this paper, we focus on a sample of UV-selected 
AGNs at z ~ 2 — 3 drawn from the Lyman break 
galaxy (LBG) (z ~ 3) an d the BX/BM samples (z ~ 2) 
(jSteidel et al1l200l l200l . The objects in our AGN sam- 
ple were identified as Type II AGNs based on narrow 
emission lines in the rest-frame ultraviolet. In such ob- 
jects, emission from the central source is obscured from 
view as the observer's sightline passes through a cen- 
tral dusty medium. This obscuration provides a better 
view of the host galaxy surrounding the central AGN. 
An analysis of the clustering of the UV-selected AGNs 
indicates that they appear to be hosted by the same dark 
matter halos as those of the parent population of non - 
AGN LBGs (ISteidel et al.ll2002l : lAdelberger et al.ll20r)5ft . 
The non-AGN LBGs, therefore, provide an ideal con- 
trol sample with which to compare the AGN LBG host 
galaxies. Here, we use stellar population synthesis (SPS) 
modeling to examine the UV through IR SEDs for both 
UV-selected AGNs and non-AGNs. SPS modeling uses 
stellar population templates to match the broad spec- 
tral energy distribution (SED) of a galaxy and therefore 
to estimate the properties of the stellar population such 
as stellar mass, age, extinction, and star-formation rate. 
There has been much effort to use SPS modeling to char- 
acterize the properties of s tar-forming ga laxies across a 
wide range in cosmic time (Shaplcy 2011L and references 
therein), and recently the method has been extended to 
model AGN host galaxies. However, using this process 
to examine AGN host properties is made difficult by non- 
stellar emission from the central source, which has to be 
carefully modeled alongside the stellar population. Host- 
galaxy SEDs have been decomposed into stellar and AGN 
contributions f or both X-ray selected Type II AGNs at 
z ~ 0.5 - 2.5 (|Mainieri et alj|2011t ISantini et all 12012ft 
and_ PR-selected AGNs at z ~ 0.7 - 3.0 (|Hickox et al.1 
12007ft . We follow a similar methodology here to charac- 
terize the UV-selected AGN host galaxy demographics 
at z ~ 2 — 3. Acc ordin gly, we build on ea rlier work from 
lErb et al.l (|2006bl ) and iKriek etail (|2007ft using a larger 
sample of AGNs and a more systematic treatment of 
AGN host galaxy stellar population modeling and selec- 
tion effects. We compare the distributions of stellar age, 
dust extinction, SFR, and stellar mass with both a non- 
AGN sample as well as a non-AGN subsample matched 
in stellar mass. We also investigate such trends as SFR 
vs. stellar mass, and U — V color vs. stellar mass, in 
order to see how the AGNs are distributed relative to 
the non-AGNs. 

One unique aspect of the UV-selected AGN sample is 
the set of existing rest-frame UV spectra for both the 
AGNs and the parent sample of non-AGNs from which 
they were selected. We can use these spectra to an- 
alyze the relationship between the host galaxy stellar 
po pulation and observed UV spectroscopic properties. 
In lHainline et al.l (|2011aft . we created a composite rest- 
frame UV spectrum for the UV-selected AGN sample, 
and studied how the composite varied as a function of 
UV spectroscopic properties such as emission line equiv- 
alent width (EW), UV spectral slope, and galaxy UV 
luminosity. Here, we use the results from SPS model- 
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ing and extend this analysis to examine how the host 
galaxy properties for this sample of AGNs are reflected 
in the UV spectra. The UV spectra can also be used to 
estimate potential selection effects as a function of host 
galaxy parameters, in order to understand the complete- 
ness of our AGN sample, and the relationship between 
black hole and host galaxy properties. 

We describe the AGN sample selection and the data 
used in the SPS modeling in $2] In Sj3j we detail our 
dual-component ("SPS+AGN") modeling technique to 
disentangle the emission from the stellar population and 
the AGN. Specifically, we use this technique on a subset 
of AGNs with IRAC data to infer how the presence of 
AGN emission affects the host galaxy parameters derived 
from SPS fitting. This analysis leads to the calculation 
of "correction factors," which are then used to model our 
full set of AGNs and robustly estimate their host galaxy 
parameters, as detailed in SJU In S}5J we compare the 
stellar properties for the AGNs to those of non-active 
star-forming galaxies, including a subsample of galax- 
ies matched in stellar mass. In detail, we compare the 
AGNs to non-active galaxies using the SFR vs. stellar 
mass relation, the rest-frame U — V color vs. stellar 
mass relation, and UV spectroscopic properties. Finally, 
in fjBJ we consider the origin and implications for the UV 
spectroscopic trends among AGNs, in terms of sample 
completeness and the connection between black hole and 
galaxy mass at z > 2. We list our conclusions in SJ7J 
Throughout this paper, we adopt the following cosmo- 
logical parameters: H n = 70 km s _1 Mpc -1 ; f2 fot , fi m , 
fl A = 1,0.3,0.7. 

2. THE UV-SELECTED AGN SAMPLE AND 
OBSERVATIONS 

The narrow-lined AGNs we model in this paper were 
selected from a parent sample of 3,059 UV-selected star- 
forming galaxies at z ~ 1.5 — 3.5. LBGs at z ~ 3 are 
chosen due to the absorption of light shortward of the 
Lyman break at 912 A by the intervening intergalactic 
medium (IGM). As this break is redshifted to optical 
wavelengths for z ~ 3 galaxies, such objects can be se- 
lected based o n their position in a U — G vs. G — TZ color- 
color diagram ()Steidel et alJl2003T ). Similar color criteria 
have been used to identify galaxies at z ~ 2, as the UGTZ 
colors of star-forming galaxies at this redshift correspond 
to a flat part of th e rest-frame UV spectrum redward of 
the L yman break (jSteidel et al.l 120041 : lAdelberge r et al.l 
120041) . The full sample of z ~ 2 - 3 UV-selected star- 
form ing galaxies is descri bed in lSteidel et al. (2003. 120041) 
and iReddv et al.l (|2008l ). From this sample, galaxies 
were identified as having an active nucleus on the ba- 
sis of features in their rest-frame UV spectrum: strong 
Lya emission accompanied by detectable emission in at 
least one other high-ionization AGN emission line, in- 
cluding N VA1240, Si IVAA 1393,1402, C IVA1549, or 
He IIA1640. These objects were selected to be "narrow- 
line," or "Type II" AGNs, as the FWHM for any of the 
emission features was less than 2000 km s _1 . The aver- 
age redshift of the 33 AGNs is (z) = 2.54 ± 0.34, with an 
average 1Z magnitude of (7Z) = 24.2 ± 0.7, and a range 
of 22.55 < 1Z< 25.72. We show the redshift distribution 
for the UV-selected AGNs in Figure [TJ We have also 
assembled a sample of UV-selected non-active galaxies, 
which have been spectroscopically confirmed at z ~ 2 — 3 
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Fig. 1. — Redshift distribution for the UV-selected AGNs (red), 
compared to simil arly-selected non-active galaxies (grey) from 
IReddv et al. (20l3). The non-AGN comparison sample spans the 
same redshift range as the AGN sample. 



l|Steidel et al.l 120031 [2001 IReddv et al.ll20Tl and have 
near- and/or mid-IR photometric coverage. This non- 
active comparison sample contains 1727 galaxies when 
limited to those in the same redshift range as our AGN 
sample (1.98 < z < 3.43) and is plotted with a grey his- 
togram in FigureQ] The non-AGN sample has an average 
redshift (z) = 2.56 ± 0.37, as compared to the average 
redshift for the AGN sample, (z) = 2.54 ± 0.34. In sub- 
sequent sections, we will use this non-AGN sample to 
highlight the impact of the AGNs on stellar population 
modeling. More details about the AGN sampl e and its 
selection are provided in Mainline et al.l (|2011al ). 

The UV, optical, and infrared photometry for this 
AGN sample comes from multiple imaging programs. 
Objects were initially selected as hi gh-redshift galaxies 
from their U, G, and 1Z magnitudes (St eidel et al.ll2003l 
2004) . The UV and optical imaging for the fields contain - 
ing our AGNs is detailed in full in IReddv etHI (|2008D . 
The galaxies that comprise this sample were selected 
from imaging characterized by a typical 3cr photomet- 
ric depth of 1Z ~ 27.0. In order to fully model these 
objects, we limited our sample to objects with near- 
and/or mid-infrared photometry At the rcdshifts we are 
studying, important spectral features that indicate age 
are shifted to the near-IR, such as the 4000 A break 
due to metal absorption lines in late type stars, and 
the Balmer break at 3646 A from A-type stars. Data 
at these and longer wavelengths are required in order 
to constrain stellar population models an d give accurate 
and robust stellar mass and age estimates (jShaplev et al.l 
l200lUReddv et al.ll2012D . For 15 of our 33 objects (in the 
Q0100, Q0142, Q0933, Q1217, GOODS, Q1422, Q1623, 
Q1700, Q2343, and Q2346 fields), near-IR J and K 
magnitudes were obtained using the WIRC instrument 
on the Palomar 200" telescope, with a typical 2" aper- 
ture 3a photometric depth of 22.6 (if), and 24.1 (J), in 
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Vega magnitudes (jErb et al.l l2006bh Fl The WIRC data 
were reduced using custom IDL scripts as described in 
IShaplev et all l|2005h . When necessary, some of the K- 
band images were smoothed to match the seeing in the 
72.-band images. Further if-band imaging for 11 ob- 
jects (in the Q0000, Q0201, Q0256, Q2233, DSF2237a, 
DSF2237b, Q2346 fields) was obtained using the PANIC 
instrument on the Magellan telescope, as described be- 
low. F or 16 objects in o ur sample, we assembled Spitzer 
IRAC (|Fazio et alll2004[ ) data in Channels 1 (3.6 /mi), 2 
(4.5 fxm), 3 (5.8 /mi), and 4 (8.0 /mi). IRAC data for our 
ob jects come from fiv e Spitzer Programs, as described 
m iReddv et all p0l2) in more detail. These data were 
reduced with custom IDL scripts that applied artifact 
correction and flat fielding, after which individual images 
were mosaiced, and drizzled. Photometry was performed 
by fitting point-spread functions (PSFs) to the images 
with positions pre-determined using optical and near-IR 
data. The full details of the I RAC reduction PSF fi tting, 
and photometry are found in IReddv et "ail (|2006bl) . We 
also assembled Spitzer /MIPS 24 /im data for 8 of our ob- 
jects, as deep surveys were undertaken in the GOODS-N 
(PI: M. Dickinson) and Westphal (PI: G. Fazio) fields, 
and the Q1623, Q1700, and Q2343 fields were imaged 
under GO 1 and 3 Spitzer programs. The MIPS data 
were flat fielded with custom IDL scripts, combined, and 
p hotometry wa s obtained wi th PSF fitting, as described 
in IReddv etUI (|2006bL l20Toh . 

As only 15 of the narrow-line AGNs in the UV-selected 
sample had existing coverage at near-IR wavelengths, 
we imaged 11 additiona l obje cts using the PANIC in- 
strument ((Martini et all 120041) on the 6.5m Magellan 
I (Baade) telescope at Las Campanas Observatory in 
September of 2008. These objects were observed over 
the course of two nights under excellent conditions (me- 
dian seeing of 0'.'5), with one additional night lost to 
weather. The data were reduced with IRAF PANIC 
reduction scripts ("p_reduce") (Martini et al. private 
communication). Each science exposure was linearized, 
dark subtracted and flat fielded, and sky subtracted us- 
ing dithered frames adjacent in time. The final images 
were distortion corrected and stacked. For the purposes 
of flux calibration, we also obtained and reduced images 
of iPersson etaLI 0998) near-IR standard stars. The typ- 
ical 3(7 photometric depths of the final images ranged be- 
tween 19.5 and 21.5 in K (Vega magnitude), with a me- 
dian of 21.0. The PANIC _ftT-band images were smoothed 
to match the seeing in the 7£-band images, which were 
typically characterized by a larger point spread func- 
tion. Phot ometry for these data w as calculated using 
SExtr actor (|Bertin fc ArnoutsllT996l ). We used SExtrac- 
tor "isophotal" magnitudes, where the isophotal area 
adopted to calculate flux in the 7?.-band image was ap- 
plied to the smoothed K-b&nd image. With this proce- 
dure, we calculated 1Z—K colors within a fixed isophotal 
area, just as the existing optical and near-IR colors had 
been estimated. 

The optical and near-IR photometric errors for all of 
our objects were estimated using a Monte Carlo ap- 
proach, where fake galaxies across a range of magnitudes 
were added to the data and recovered using SExtractor, 
with the same input parameters as those used to calcu- 

3 24.4 (K) and 25.0 (J) in AB magnitudes 
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Fig. 2. — Fractional 1Z—K color distribution for our sample of 
z ~ 2—3 AGNs, compared w ith similarly-selected non-active galax- 
ies from IReddv et al.l pOT2f). The 11- K colors for AGNs are, on 
average, redder than for non-active galaxies, indicating older stellar 
populations and/or more dust extinction. 

late the actual photometry. This method resulted in dis- 
tributions of the differences in input and recovered mag- 
nitudes in various bins of recovered magnitude. The error 
that we assigned to each photometric measurement was 
the standard deviation in the relevant bin. Photometric 
errors for the IRAC data were obtained with a similar 
method, where the photometry was calcul ated using a 
custom PSF fitting method as described in IReddv et all 
(2012). The full UGKJK+ Spitzer photometric datasets 
for our objects are listed in Tabic Q] 

3. SED MODELING 

Before discussing the stellar populations of our AGN 
hosts, we examine their observed colors. The measured 
1Z—K distribution for our full sample of AGNs is shown 
in Figure [2] along with the 1Z—K c olors for non-activ e 
t/GTvL-selected star-forming galaxies (jReddy et al.ll2012T) . 
Although these AGNs and non-active galaxies were se- 
lected using the same rest-UV photometric criteria, the 
redder 1Z—K colors for the AGNs may be indicative of a 
difference in the host stellar populations of active and 
non-active galaxies. Specifically, these colors imply a 
more evolved stellar population and/or higher dust ex- 
tinction for the AGN host galaxies, and we can use SPS 
modeling to interpret the color information for our sam- 
ple of active galaxies in terms of physical properties. 

In SPS modeling, the observed UV through infrared 
photometry of a galaxy is fit by assuming a population 
of stars governed by a set of parameters, including star- 
formation history, metallicity, initial mass function, stel- 
lar age and mass, and degree of dust extinction. Stellar 
population models use a stellar library of the observed 
spectra of stars of different spectral types, as well as a de- 
scription of the evolution of stars of different masses and 
metallicities to produce integrated spectra for a popula- 
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tion of co-evolving stars at various ages and for a given 
star-formation history. These spectra are passed through 
photometric filters to predict galaxy SED luminosities 
and colors, which are compared to observations to deter- 
mine the physical properties of the host galaxies. 

The SPS approach to SED fitting has been widely 
used to describe the stellar popula ti ons of high- 
redshift galaxies dSawicki fc Yed 19981: iShaplev et al.l 
l200lt iPapovich et al.ll2001l ; lShaplevll2011l and references 
therein). For our sample of active galaxies, however, we 
must use caution because of the effect that non-stellar 
continuum and line emission may have on the fitted pa- 
rameters. In order to constrain the effect of AGN emis- 
sion on the stellar population synthesis modeling pro- 
cess, we modeled the observed SED with both a stellar 
population component as well as an AGN component 
("SPS+AGN"). In this approach, the observed colors 
and magnitudes are decomposed at each wavelength into 
a sum of the fluxes from the underlying stellar population 
and emission excited by the central AGN. Our objects are 
narrow-lined AGNs, and thus the non-stellar emission 
should be limit ed to the rest-fra me near- to mid-infrared 
wavelengths (|Assef et al.l 120101 ). At the same time, we 
can use our dual-component modeling for objects with 
photometric information at infrared wavelengths both 
to understand the true stellar populations of AGN host 
galaxies, and also to understand how the presence of an 
AGN affects the best-fit stellar parameters. 

In this section, we describe our dual-component model- 
ing in detail. We start by discussing our sample of AGNs 
with IRAC data, which we can use to describe the AGN 
contribution to the SED (i j3.ll) . We then use our dual- 
component modeling procedure to calculate corrections 
to host galaxy parameters that account for broadband 
AGN emission ( fc|3.2[) . We also correct our photometry 
for the presence of strong, discrete AGN emission lines 
( M3.3p . and conclude with a section detailing our system- 
atic uncertainties 



3.1. IRAC Data 

As described in S}3J there is UGTZJK photometry for 
the majority of our sample, but only 16 AGNs have IRAC 
measurements. Of these 16 objects, 11 are characterized 
by "power-law" emission in the rest-frame near-IR, where 
the fluxes in the IRAC bands increase monotonically 
towards longer wavelengths. This emission is hypothe- 
sized to arise in an AGN from therma l or non-thermal 
emission near the dusty active nucleus dRees et al.lll969t 
iNeueebauer et al.l 119791: lElvis et al.l I1994I ). An addi- 
tional two objects with IRAC coverage (Q1623-BX454 
and Westphal-MM47, only including Channels 1 and 2) 
have flat mid-IR slopes, indicating a weak AGN as com- 
pared to the stellar population. As these objects do not 
have if-band data, we fit their SEDs solely with a stel- 
lar population. The final three objects (Q1217-BX46, 
Q1623-BX747, and Q2233-MD21) only have IRAC data 
in Channel 1 , which was not sufficient for the purpose of 
discerning the presence of a power-law in the rest-frame 
near-IR. We fit the SEDs for these objects, including the 
IRAC Channel 1 data, with a stellar population only. 

To place our sample of UV-selected IRAC AGNs in 
context, we compare it to other galaxy samples with cov- 
erage in all four IRAC channels. In Figure [3j we plot the 
11 objects with IRAC power-law slopes on the color-color 



diagrams from lLacv et al.l ()2004D and lStern et al.l (|2005|) . 
which were initially used to select galaxies with red IR 
colors as possible AGNs. The AGNs with detections in 
all four IRAC bands are plotted as red points, and, in the 
left figure, we use dashed lines to indicate the positions of 
three AGNs with only Channel 4 and Channel 2 photom- 
etry. In this left panel , we show the AGN selection box 
from lLacv et al.l (|2004D with long d ashed lines. Each o f 
our AGNs lies in the detection area. iDonlev et all ([20121 ) 
revised the selection criteria to account for contamina- 
tion due to star formation, and we indicate the revised 
AGN selection area with a dark line. All but one of our 
objects fall within this revised area. On the right, we 
show the AGN selection criteria from lStern et al.l (2005*) 
with long dashed lines, and all but one object falls within 
this area. HDF-BM Z1156 falls outside of the lStern et al.1 
and IDonlev et al.1 selection areas (alt hough it is co nsis- 
tent with being included within the IDonlev et al.l se- 
lection area when its photometric uncertainty is taken 
into accountj^- The grey points on this diagram are 
those objects from the Spitzer First Look Survey with 
clean detections in all four IRAC Channels across all red- 
shifts, and illustrate the full range in color-colo r spac e 
spanned by extragalactic sources (jFadda et al.l [2004). 
We have also plotted those objects from our non-active 
UV-selected z ~ 2— 3 comparison sample with IRAC cov- 
erage in all four channels. These points are shown as dark 
grey squares. Quite strikingly, UV-selected non-AGNs 
and AGNs occupy largely distinct regions in IRAC color- 
color space. Only a small fraction of the UV-selected 
non-AGNs fall within the Donley et al. (2012) selection 
area, and these objects may comprise a sample of ad- 
ditional AGNs without strong high-ionization emission 
lines in their rest-frame UV spectra. 

To compare to high-redshift power-law galaxies, we 
have plotted a sa mple of z > 1.4 power-law sources from 
iPark et all (|2010D on both diagrams as blue points. Our 
active sample spans a large range in IRAC power laws. 
The specific values for the power-law slope for our objects 
(a, where f u oc i> a ) range between —0.62 and —3.50, with 
an average (median) of a = —1.75 (—1.79). These val- 
ues sh ow the same range as those presented in lPark et al.1 
(2010), although the objects in our sample with the most 
negative slope values, Q1623-BX663 and Q1700-MD174, 
would be among the reddest of the Spitzer power-law 
galaxies (the power law for Q1700-MD174 is fit from 
only the IRAC Channel 2 and Channel 4 points). Our 
galaxies are thus representative of the full range of high- 
redshift Spitzer power-law galaxies, and should illustrate 
AGN properties for a range of SED types. It is notable 
that the majority of the UV-selected AGNs in our sam- 
ple with IRAC data have power-laws, supporting their 
id entification as AGNs , consistent with the observation 
in IReddv et all lf2006a) that all but one UV-identified 
AGN in the GOODS-N field had significant excess flux 
at 8 /im. 

3.2. Stellar Population and AGN Modeling Procedure 

4 Wc note that HDF-BMZ1156, fo rmerly known as H DF- 
oMD49, has previou sly been discussed in Steidel et al. ( 20021) and 
IReddv et at] J2006ar) . This object was only weakly detected in the 
Chandra 2 Msec data, an d has been identified as an example of 
a Compton-thick AGN by Alexander et al. (2008b), based on its 
combined mid-IR and X-ray properties. 
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but HDF-BMZ1156 would be selected as AGNs u nder the IStern etaH 1 12003 ) and lDonlev et al.l ( 120121) criteria. Predictions of AGN colors 
as a function of redshift from Donlcv ct al. (2012) find that strong AGNs inhabit the space occupied by HDF-BMZ1156 at z ~ 2 — 3. 



In order to accurately describe the SEDs for our sam- 
ple, we require a suite of stellar population models and 
an AGN template. For our sampl e of z ~ 2 — 3 AGNs 
we chose iBruzual fc Chariot] (|2003l . hereafter BC03) stel- 
lar population models. Recently, much attention has 
been focused on the degree to which thermally pulsat- 
ing asymptotic giant branch (TP-AGB) stars contribute 
to the near-IR luminosity of a population of . 5 — 2 
Gyr age stars (iMarastonl l2005t iMaraston etaH 120061 : 
lEminian et"aDl2008f ). Different treatments of TP-AGB 
stars systematic ally affect derived stellar mass and ages 
(|Marastonll2005l ) . In order to gauge how the presence of 
TP-AGB stars affects our inferred stellar para meters, we 
also a pplied the stellar population models of Maraston 
(|2005l hereafter Mar05) and Chariot & Bruzual (2012) 
(private communication, hereafter CB12) to our SEDs. 
For the BC03 and CB12 models, we adopted the Padova 
1994 stellar evolut ion track s. We s elected solar metallic- 
ity models and a iChabrierl (|2003f ) initial mass function 
(IMF) extending from 0.1 to 100 M . We adopted so- 
lar metallicity models as non-active galaxies in a simi- 
lar mass range a ppear to have meta llicities only slightly 
lower than solar (jShaplev et al.1 2004) . It should be noted 
that using a Salpeter IMF increases the best-fit stellar 
masses and star-formation rates by a factor of 1.8 with 
little effect on the predicted spectral shape. To ac count 
for dust extinction, we used a lCalzetti et al. (2000) star- 
burst attenuation law, which on average provides an ac- 



curate description of the reddening and attenuation of 
the UV stellar con tinuum in both nearby and distant 
starbu rst galaxies (jReddv fc Steidell 12004 iReddv et"ai] 
I2006bl ). 

For the BC03 and CB12 models, we use constant star 
formation (CSF) models. For the Mar05 models, we 
used an exponentially declining SF model of the form 
SFR(i) oc exp(— t/r), with an e-folding time of r = 20 
Gyr as the closest approximation of a CSF model. We 
chose to not model with more complex star-formation 
histories, such as exponentially declining star-forming 
models, short-duration bursts superposed on top of an 
exponentially declining star-formation history, or expo- 
nentially incr easing star fo r matio n models. Recent work 
presented in IReddv et al.l (|2012l ) indicates that expo- 
nentially declining star-formation models result in SED- 
derived SFR values that are lower on average than those 
derived from a combination of UV and IR dat£0- Ex- 
ponentially increasing star formation models have been 
shown to apply to galaxies at z > 2, and naturally 
arise in hydrodynamic and semi-analytic models of star 



5 Reddy ct al. derive the SFR/^j+jyv' by combining SFRs 
measured from photometry probing the rest-frame UV and IR. 
Observed-frame optical photometry is used to calculate the rest- 
frame UV (1700 A) luminosity, and MIPS 24 /im photometry is 
used to calculate the rest-frame IR (8 fim) luminosity. These lu- 
minosities are converted to SFRs assuming typical galaxy ages and 
star-formation histories, and then combined to form SFR//j_|_[/^. 
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Fig. 4. — Histogram of reduced \ 2 values from SPS-only model 
fits to UV-selected IRAC AGNs (red), and UV-selected non-AGNs 
(grey). The IRAC AGNs are not well- fit by SPS-only models due 
to emission in the IRAC Channels, while the bulk of the non-AGNs 
are well fit ((xfLd) = ^-3 i 1.8). In the inset, we show the best-fit 
SPS-only SED model for an example AGN, HDF-MMD12, using 
a BC03 CSF model. The power law seen in the IRAC bands is 
indicative of the presence of an AGN, which results in a poor fit 
(Xred = for a stellar emission-only model. The errors on the 
IRAC points are smaller than the size of the data points. The solid 
curve is the fit derived taking into account all of the data. 

formation, and are suggested by the relationship be- 
tween SFR and stellar mass observed at high redshift. 
(IMaraston et al.ll20ia|Papovich et aT1l2011l:lRe"ddv et al.l 
2012?). These histories give SFRs that are similar to the 
SFRs derived from CSF models, along with older ages 
HReddv et al.l[20L2T) . 

Stellar population models, even those with large con- 
tributions from old, red stars, cannot successfully model 
the infrared emission seen in an IRAC power-law object. 
In Figure 21 we show in red the distribution in reduced 
X 2 (xled) f° un d by modeling our IRAC AGNs with only 
a stellar population. We also plot the xt e d distribution 
for the non-AGN comparison sample in grey, using the 
subsample of 231 non-AGNs with SED coverage extend- 
ing out to 8/im. The AGNs have higher Xred values than 
the bulk of the non-AGN comparison population, indi- 
cating poor fits. For the non-AGNs, (xled) = 2 - 3 ± L8 > 
which is lower than all but 1 of our IRAC AGNs (Q0142- 
BX256, which does not have IRAC Channel 4 data, but 
has large errors on the Channel 3 data point). In the 
inset we show a CSF fit for one of our objects, HDF- 
MMD12, an AGN at z = 2.648. The best fit is very poor, 
with a xfed — 36, due to the presence of the strong AGN 
power-law emission in the IRAC bands. Often, the SEDs 
of Type II AGN hosts are modeled excluding points red- 
wards of the rest-frame optical, as the AGN emission for 
narrow-lined objects sho uld be minimal at rest-frame UV 
and o ptical wavelengths (B arger et al . 2005; B undv et al.l 
120081) . With our sample of IRAC AGNs, we can model 
the power-law AGN emission and the stellar emission si- 



multaneously to see exactly how the presence of an AGN 
affects the best-fit stellar population parameters. We 
also note that there exists a tail of galaxies without signs 
of AGN activity in their rest-frame UV spectra with large 
Xred va l ues from SPS-only model fits. We will discuss 
these objects in f]4J] and in a forthcoming paper. 

We then require a prescription for the non-stellar AGN 
emission in our dual-component modeling. Often, the 
non-stellar continuum contribution to the near-IR emis- 
sion is removed by assuming that this portion of the 
spectrum can be represented by a power law of the 
form f v oc y~ a , where f v is the flux density at a 
freque ncy v (|Alonso-Herrero et al.l 120031 : lHainline et al.l 
1201 lbl) . Since this method of modeling the AGN emis- 
sion has no explicit basis in the phys ics of the AGN, we 
instead opt for the AGN model from lAssef et all (|2010f) 
to describe our AGN infrared emission^. The Ass ef et all 
AGN model was created using an it erative method, start- 
ing with the IRichards et al.l (|2006l ) average quasar SED 
template, modified to better replicate the behavior in the 
UV predicted by an accretion disk model (in an attempt 
to minimize, in part, contamination by host galaxy emis- 
sion). This template was applied to the SEDs of AGNs 
from the NOAO Deep Wide-Field Survey in the Bootes 
field, and iteratively updated to fit the full sample of 
data. The final AGN template is thus physically moti- 
vated, and, as the template has been both normalized 
to a distance of 10 pc and created to have an integrated 
luminosity of 10 10 Lq, its application to our data gives 
us a n estimate o f the b olometric luminosity of the AGN. 
The lAssef et all (|2010D AGN template describes a Type 
I AGN, with strong emission in the rest-frame ultraviolet 
thought to be emitted from the accretion disk, as well as 
emission in the rest-frame near- to mid-IR from warm 
and hot dust emission. At the sam e time, the AGNs 
that we model with the lAssef et al.l template were se- 
lected to be Type II by virtue of narrow emission lines 
seen in the rest-frame UV, where accretion disk emis- 
sion should be obscured by dust. To model our Ty pe II 
AGNs, following the analysis of lAssef et al.l (|2010l) . we 
apply dust extinction to the AGN template using a red- 
dening curve derived from the Small Magellanic Cloud 
for A < 3300A ( (Gordon fc Clavton ||1998D and following 
the G alactic curve at longer wavelengths (Cardclli et al. 
1989). The resulting template should then consist of 
strong IR emission from hot dust near the active nucleus, 
which we use to model the IRAC power-law. 

Since the objects we are modeling were selected to 
be Type II AGNs, we made a prior assumption of 
E{B — V)agn > L0. We model these objects by taking 
a Type I AGN template and applying extinction in or- 
der to produce a Type II AGN template to fit the SEDs. 
For the Type II AGNs discussed here, the rest-frame UV 
and optical light should be dominated by stellar emis- 
sion, given the narrow emission lines in the rest-frame 
UV s uggesting signif i cant o bscuration of the central en- 
gine (jHainline et al.l 12011a!) . Low values of AGN ex- 
tinction, on the other hand, would result in SEDs with 



6 We also explored t he use of multiple diff erent empirical Type 
II AGN t emplates fro m lPolletta et all (1200711 . but found that the 
reddened Asse fet al.l model was better able to reproduce the 
observed phot o metry of our UV-selected AGN sample (see also 
Mainicri ct al. (2011) for a similar methodology). 
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prominent AGN continuum emission even at rest-frame 
UV wavelengths, as well as broad-line region emission. 
Although the initial modeling results for two systems, 
HDF-BMZ1384 and Q2343-BX333, indicated best fit val- 
ues of E{B — V)agn less than 1.0, we rejected such fits 
on the basis that these are Type II objects. Furthermore, 
the x 2 values for these low E(B— V)agn fits were not sig- 
nificantly lower than for fits with E{B — V)agn at values 
similar to the bulk of the sample (E(B — V)agn ~ 2 — 7). 
T he same lower l i mit o n E(B — V)agn is also adopted 
in iMainieri etHl (|2011[ ) for a similar SPS+AGN fitting 
procedure used to fit X-ray selected Type II quasars 
at z > 0.8. Throughout this paper, we report only 
those fits where E(B — V)agn > 1-0. We also limit 
E(B - V)agn < 8.0, as larger values for the AGN ex- 
tinction did not improve the fits for the objects we ex- 
amined. 

We used the lAssef et al.l AGN template, along with 
the stellar population models, to fit the observed SEDs 
of our IRAC subsam ple. The overa ll procedure we use d 
is similar to that of iShaplev etail (|200l I200H I2001 . 
but here the predicted colors are the sum of stellar emis- 
sion and non-stellar emission from our AGN model. The 
best-fit model therefore yields constraints on both stel- 
lar population as well as AGN parameters. The stel- 
lar parameters that we modeled included dust extinction 
[parameterized as E(B — V)sf], stellar population age 
{tgp), star-formation rate (SFR), stellar mass (M t ), and 
in the case of the Mar05 models, star-formation history 
(r). The ages that we used ranged between 50 Myr and 
the age of the universe at the redshift of the modeled 
galaxy. This lower limit of 50 Myr was chosen to exclude 
models with ages younger than the dyn amical timescales 
estimated for ty pical z ~ 2 - 3 LBGs (|Erb et alj|2006at 
ILaw et all 120071 ). We considered stellar extinctions be- 
tween E(B-V) S F = 0.0 and CiQ. The AGN parameters 
included a central source extinction, E(B — V)agn, and 
a normalization constant, Aagn- 

In the first step of the SED fitting procedure, we cre- 
ated total SPS+AGN models for each potential combina- 
tion of input parameters. The stellar model flux at each 
wavelength (with a specific star-formation history, age, 
and extinction value) was summed with the AGN model 
flux (with a specific AGN extinction value, and multi- 
plied by the normalization factor Aagn)- In order to add 
the AGN template to our SPS models in a meaningful 
way, we used values of Aaon that reflect ed how the SPS 
models are normalized. The lAssef et al.l AGN template 
is normalized to an AGN bolometric luminosity L{, i = 
1O 1O L0, where this luminosity is calculated by integrat- 
ing the entire template longward of Lya. In terms of 
the SPS models, the BC03 and CB12 CSF templates are 
normalized at 1 M Q yr _1 of star formation. Thus, for the 
BC03 and CB12 CSF fits, N AG n = L boi /{SFR x 10 10 ). 
The Mar05 models are all normalized to 1 M Q of stars, 

7 It is important to note that SPS-only models with arbitrarily 
high values of E(B — V)sf do not accurately fit the UV-selected 
AGN SEDs. Due to the specific shape of the SED of a star-forming 
stellar population, the application of high values of E(B — V)sf 
cannot not simultaneously replicate the rest-frame UV and optical 
photometry as well as the rest-frame near-IR power-law that is 
typical of our UV-selected AGNs. In the best-case scenario, such 
fits vastly underpredict the data at shorter wavelengths while only 
roughly fitting those in the rest-frame near-IR. 



and so Nagn = L bo i/(M* x 10 10 ) for these fits. While we 
report the best-fitting Nagn values, it is useful also to 
compare L bo i between fits, which indicates the strength 
of the AGN for each object. 

The combined AGN and stellar model spectrum was 
further attenuated b y IGM absorption from neutral hy- 
drogen (Madau 1995), which only affects the predicted U 
and/or G magnitudes, depending on the redshift. These 
output SPS+AGN templates for a given set of param- 
eters were compared in magnitude space to the optical 
through mid-IR SED for a specific object (though, not 
including Spitzer MIPS data, see ^4.1[) . and a value of x 2 
for the fit was calculated. In these fits, we only modeled 
photometric detections. In a few cases, there were upper 
limits in IRAC Channels that were not modeled (Channel 
4 for Q0142-BX256, Channel 3 and Channel 4 for Q1623- 
BX454, Channel 2 for Q1623-BX74, and Channel 3 for 
Q2233-MD21). In all cases except that of Q0142-BX256, 
the IRAC upper limits are entirely consistent with the 
best-fit SED model. We sought to scale our total mod- 
els to the observed SEDs such that the calculated \ 2 was 
minimized, and this normalization naturally led to an es- 
timation of the star-formation rate and stellar mass. The 
parameters [E(B-V)sf, tSF, SFR, M*, E(B-V) A gn, 
and JVagn] resulting in the minimum \ 2 with respect to 
the measured photometry for an object are referred to as 
the "best-fit" parameters. 

3.3. Correcting the Photometry for the Presence of 
Strong Emission Lines 

Flux contribution from strong nebular emission 
lines excited by the presence of an AGN can be 
a signi ficant source of con t amina tion in broadband 
SEDs dSchaerer fc de Barrosl 1201 lh . As described in 
lHainline et al.l (|2011a|) . the rest-frame UV spectra for the 
z ~ 2 — 3 Type II AGNs described here include emis- 
sion lines such as Lya, with an average rest-frame EW 
of 80 A. This strong emission line, as well as others in 
the rest-frame UV, such as C IVA1549 and He IIA1640, 
can enhance the flux in the U or G bands (depend- 
ing on redshift), and alter the observed U — G or 
G—1Z colors. Rest-frame optical features such as Ha or 
[O III]A5007 can similarly add flux in the near-IR, and 
subsequently affect the measured 1Z—K colors. We adopt 
a method of correcting for these emission lines outlined 
in lPapovich et al.l ()200lD . where a rest- frame EW, EWo, 
for an object at a redshift, z, introduces a change in the 
magnitude: 



Am = —2.5 log 



1 



EWo(l + z) 



AX 



(1) 



In this equation, AA is the width of the filter bandpass. 
We calculated Am for each measured emission line, and, 
where its absolute value was larger than the photomet- 
ric error, we corrected for this magnitude difference to 
produce the final photometry. 

Strong emission lines do not affect the U or G bands 
for 19 out of 33 AGNs in our sample at 2.17 < z < 2.48. 
For the remainder, w e used UV spectra (described in 
Hainl ine et al.l l2011a|) to measure EW values for the 
strong emission lines such as Lya, C IV, and He II. Lya 
was the chief contaminant to the G band for this redshift 
range, with the average (median) value of the magnitude 
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difference Am = 0.36(0.34). The other UV lines were 
not strong enough to significantly contaminate the pho- 
tometry: in objects where C IV was detected, the EW 
was on average (median) 48% (17%) of the EW of Lya, 
and similarly, in objects where He II was detected, the 
EW was on average (median) 25% (13%) of the EW of 
Lya. However, for one object (HDF-BMZ1156), very 
strong C IV and He II contaminated the G band (at the 
rcdshift of this object, Lya fell bluewards of the G band), 
with a combined Am = 0.37. 

The K band contains Ha for 20 out of 33 AGNs, and 
[O III]A5007 for an additional 4 objects. Near-IR spectra 
have only been obtained for a subset of our AGNs, as 
part of an ongoing project to measure the rest -frame op- 
tical s pectroscopic properties for this sample (jErb et al.l 
2006b, Hainline et al., in prep). For this analysis, there 
are 9 objects with near-IR spectral coverage (8 objects 
with Ha measurements, and one with a measurement of 
[O III] ) , and we have determined the Ha and [O III] fluxes 
for the available lines. These spectra do not show a sig- 
nificant continuum, so we used the if -band magnitude in 
order to estimate the continuum flux density (corrected 
for line flux), and therefore the EW for Ha and [O III]. 
The average (median) magnitude difference from the Ha 
emission lines is Am = 0.18(0.12). It should be cau- 
tioned that the majority of our objects do not have rest- 
frame optical spectra, and Ha or [O III] could contribute 
to the K band photometry for these AGNs. However, 
we expect that the effect should not be large compared 
to the if -band uncertainties, based on objects with such 
measurements . 

3.4. Systematic Uncertainties 

While the uncertainties relating to single-component 
stellar population fitting are well described in 
IShaplev et all (|2005f) . the addition of an AGN compo- 
nent to the fits introduces new model degeneracies and 
systematics. For example, without an AGN component, 
the inferred extinction, E(B — V), depends on the chosen 
star-formation history and age. For CSF models, a given 
set of galaxy G ~ 1Z and 1Z—K colors can be produced 
by models with older ages and less dust extinction, or 
younger ages and more dust extinction. 

With the addition of an AGN template to the SPS 
modeling, we find covariances between AGN model ex- 
tinction {E(B — V)agn) and normalization (Nagn), 
as well as between the AGN parameters and the stel- 
lar parameters. Within the IRAC sample described in 
t j3.11 a given set of photometric values for our individ- 
ual AGNs can typically be explained by model fits with 
lower values of E(B — V)agn along with lower values 
of iVAGN, or higher values of E(B — V)agn along with 
higher values of A^agn ■ This degeneracy has an effect on 
the best-fit stellar parameters. For a given object, lower 
E{B — V)agn and ./Vagn model fits result in lower val- 
ues of E(B — V)sf and correspondingly lower SFRs. As 
described above, lower E(B — V)sf fits are also accom- 
panied by older ages. We can explain these covariances 
qualitatively. For a given intrinsic AGN power-law, lower 
values of E(B — V)agn indicate more AGN contamina- 
tion at shorter wavelengths, and these fits have low values 
of Nagn in order to fit the observed SED without over- 
predicting the UG1Z flux. The increased flux at shorter 
wavelengths from a low E(B — V)agn AGN template 



is accompanied by lower E(B — V)sf values, as the in- 
creased AGN flux requires less dust to account for the red 
UV to optical continuum (1Z—K color). Our modeling 
procedure works by adding the AGN template (modi- 
fied by some dust obscuration, and multiplied by A^agn) 
to a stellar population model and then normalizing the 
result to the observed SED to determine the SFR. As 
a low- E(B — V)agn, Iow-A^agn dual-component model 
has more flux from the AGN at shorter wavelengths, such 
a model predicted that young stars contribute less flux 
in the rest-frame UV, which results in lower SFR values. 
We demonstrate this covariance using Q0100-BX172 as 
an example in Figure [5] 

In order to account for covariance among the best-fit 
stellar and AGN parameters, errors in best-fit parame- 
ters were estimated using a Monte Carlo analysis. We 
created a large set (N = 500) of fake SEDs by varying 
the observed colors and magnitudes in a manner consis- 
tent with the photometric errors. Each fake SED was 
modeled using the same method as for the actual data, 
where the best-fit stellar and AGN parameters were those 
that resulted in a minimization of \ 2 - We estimated the 
la errors for the best-fit model parameters by calculat- 
ing the standard deviation of the distribution on each 
parameter from our Monte Carlo simulation results. 

4. STELLAR POPULATION MODELING RESULTS 
4.1. AGN and Stellar Population Modeling Results 

The dual-component SPS+AGN fits successfully re- 
produce the range in SEDs found in the IRAC AGN 
subsample. We show the SED fits for the 11 IRAC 
AGNs in Figure [6l In each panel, the grey curve rep- 
resents the best-fitting dual-component model, which is 
decomposed into the sum o f fluxes from a stellar popu- 
lation (blue curve) and an lAssef et al.l AGN template 
(red curve). Also shown are the x 2 and xled values 
output from the modeling program. The poorest fits, 
with the highest values of Xredi are f° r t nose objects 
where the if -band flux was greater than predicted from 
the fit, as in HDF-BMZ1156, HDF-BMZ1384, and to a 
lesser extent Q1700-MD174 and Q2343-BX333. The un- 
derprediction of if-band flux could be a result of the 
choice of star-formation history, or the contamination 
from Ha emission to the flux in the if-band. It should 
be noted that we have a rest-frame optical spectrum for 
HDF-BMZ1156, and have made a correction for the pres- 
ence of Ha in this object (Am = 0.07), and yet the 
best-fitting model still underpredicts the K band. We 
have also calculated the fraction of the total emission in 
the best-fit SPS+AGN model at rest-frame 1 /xm that 
is contributed by the AGN template. These values of 
fAGN,ifim are reported in the bottom right corner of each 
panel in Figure 6, and express the relative strength be- 
tween the emission from the AGN and the stellar pop- 
ulation. The values of Jagn i^m range from < 0.01 
for HDF-BX160 and Q1700-MD157 to 0.79 for Q1623- 
BX663. The best-fitting parameters for the IRAC AGNs 
are reported in Table [U The IRAC AGNs have average 
(median) values and standard deviations for the best- 
fit parameters of (log(Mass)csF) = 10.82 (10.76) ± 0.22, 
(SFRcsf) - 56 (37) ± 51 M yr" 1 , (E(B - V) S f,csf) = 
0.23 (0.24)±0.08, and (log(Age)csF) = 3.23 (3.21)±0.25. 

Based on the range in best-fitting parameters, we can 
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Fig. 5. — Confidence intervals for SED fits to Q0100-BXf72, showing the covariance between Nagn and E(B — V)agn in our fitting 
results. This plot shows the density of points from the Monte Carlo error results, where darker shading indicates a higher density of points. 
The black star is the best-fit value from the SED fitting as shown in Table |21 The majority of fits show low values for Nagn across a range 
of E(B — V)agn, with both increasing together proportionally at larger values. Below the confidence interval plot, we show two sample 
fits for the black squares, indicating how the stellar parameters (with units of log(Agc) in Myr, SFR in Mq yr - x , and log(Mass) in Mq) 
change with Nagn and E(B — V)agn- The blue line is the best-fit stellar model (BC03), while the red line is the best-fit AGN model. 
The sum of the models in each panel is given as the grey line. 



split the IRAC AGN host galaxies into two popula- 
tions based on the slope of the IRAC power law (Fig- 
ure 13]). The first population contains those objects in 
our IRAC subsample with shallower observed power-law 
slopes. Most of the galaxies that comprise this pop- 
ulation (Q2343-BX333, Q0100-BX172, Q0142-BX195, 
HDF-BX160, and Q1700-MD157) have dual-component 
fits where the AGN only becomes the dominant source 
of emission at 6 — 8/im observed- frame (Fig. |6]). These 
objects are best fit by models with low values of Nagn 
along with a range of E(B — V)agn values. The one ex- 
ception is HDF-MMD12, which has a low value for Nagn 
accompanied by strong AGN emission at observed-frame 
wavelengths shorter than 8 /mi. HDF-MMD12 is best fit 
by a relatively young stellar population (tsF ~ 200 Myr). 
This model SED has a weak Balmer break with respect 
to those of the other objects, and less flux from stars 
at lon ger wavelengths, in comparison to the Asse t et al.l 
(2010) AGN template. 
The second population of IRAC AGNs consists of 



objects with much steeper (redder) power-law slopes. 
These objects (HDF-BMZ1384, HDF-BMZ1156, Q1700- 
MD174, Q0142-BX256, and Q1623-BX663) are best fit 
by models with larger values of Nagn, in which the AGN 
emission already becomes dominant at wavelengths as 
short as 5 — 6/j.m in the observed frame. Correspond- 
ingly, the rest-frame UV AGN emission lines in the ob- 
jects with steeper power-law slopes are stronger on aver- 
age than those measured in the shallow power-law ob- 
jects. The AGN bolometric luminosities of the steep 
power-law objects are also larger, on average, than those 
of the shallow power-law objects. Q1623-BX663 is the 
most extreme example. This galaxy has a steep power- 
law (a — —2.24), with a large value of Nagn accompa- 
nied by a low value for the AGN extinction. As AGN 
emission in this object is dominant even to the K band 
for CSF models, Q1623-BX663 represents a "worst-case 
scenario" for modeling the AGNs using UG1ZJK pho- 
tometry alone and not including an AGN component. 
Q1623-BX663 has one of the steepest slopes observed 
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Fig. 6. — SED models for the 11 objects with IRAC power-law fits. The blue line is the best-fit stellar model (BC03), while the red line 
is the best-fit AGN model. The sum of the models is shown as the grey line. The best-fit parameters for these models are given in Table |21 
and the log(Mass) shown is in units of solar masses. In the bottom right corner, we show the fraction of AGN flux to total flux measured 
from the best-fitting models at rest-frame 1 /im. Values listed as < 0.01 have a negligible AGN contribution at 1/im. 

in our IRAC subsample, and it is unlikely that many 
objects in the full AGN sample will have similar AGN- 
dominated fits. Together, these AGN populations show 
a wide range of AGN contamination, and should provide 
a representative subset for correcting the full sample for 
AGN emission. 

As a sanity check, we apply our dual-component mod- 
eling procedure to the subsample of non-AGNs (dis- 
cussed in [J2] and ^3.2|) that have coverage extending to 
8 \im. Half of these objects have best-fit models with 
Nagn = 0.00, and the full distribution is highly skewed 
towards either zero or extremely low values of Nagn, 
in contrast to our UV-selected AGNs with IRAC cover- 
age, all 11 of which are fit by models with a non-zero 
value for Nagn- We assessed the probability of draw- 
ing a sample of 11 non-zero values from the non-AGN 
Nagn distribution, using 10 6 random trials. This test 
revealed that the likelihood of randomly drawing such 
a sample from the non-AGN distribution is only 0.05%, 
again demonstrating the distinct nature of the SEDs of 
the UV-selected AGNs. We do note that there is a small 
population of objects in the non-AGN sample that is best 



fit by models with Nagn significantly greater than zero. 
While a discussion of these objects is outside of the scope 
of the current work, their strong IRAC power-law emis- 
sion indicates that they may be hosting an AGN, but do 
not show significant emission lines in their UV spectra. 
These objects will be followed up in a paper analyzing 
LBG/BX/BM AGNs selected on the basis of their in- 
frared emission. As these objects comprise only a small 
fraction of the total non-AGN sample, they will not sig- 
nificantly affect our conclusions regarding the distribu- 
tion of non-AGN stellar populations. 

As discussed in £ 13 . 2 1. we also used the CB12 and Mar05 
models in our dual-component SPS+AGN fitting. Here 
we report the differences in derived parameters for the 
CSF models. For the CB12 modeling, the stellar masses 
and ages from the fitting are typically 70% of those de- 
rived using the BC03 models, while the CB12 SFR and 
E(B — V)sf values are statistically equivalent to the 
BC03 values. The Mar05 models produced similar dif- 
ferences. For the Mar05 modeling, the stellar masses 
from the fitting are typically 60% of those derived using 
the BC03 models, and the ages are 50% of those derived 
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using BC03. The E(B — V)sfr values derived using 
both models were statistically the same, and the Mar05 
SFR values were larger than the BC03 values by 25%. 
Despite these small systematic differences among stellar 
population models, our main conclusions regarding the 
AGN contamination do not depend on the specific stellar 
population model adopted. 

We have MIPS photometry for eight of our objects 
with IRAC data. We do not use the MIPS data to con- 
strain our SPS modeling, and find that for the majority 
of the objects with MIPS data, the best-fit AGN tem- 
plate overpredicts the MIPS flux in our objects by a fac- 
tor of two on average. Examining the range of stellar 
and AGN parameters that best fit our optical through 
IRAC data, we find that attempting to fit for the MIPS 
data leads to best-fit models that systematically under- 
predict the IRAC fluxes, due to the shape of the reddened 
lAssef et al.l templa te. Analysis of of the mi d-IR SED of 
HDF-BMZ1156 by lAlexander et al.l (l2008bfl sheds light 
on the origin of this discrepancy. lAlexander et all find 
that a Type II AGN template (i.e., the SED of NGC 
1068) provides a good fit to the photometry of HDF- 
BMZ1156 at wavelengths from 16 through 70 /im. How- 
ever, an additional hot (~ 1000K) dust component, of- 
ten fo und in obscured AGNs (e.g., Alonso-Hc rrero et al.l 
2001), is required to fit the IRAC datapoints at shorter 
wavelengths. If the other AGNs in our sample are sim- 
ilar to HDF-BMZ1156, our use of a single component 
to describe both IRAC and longer-wavelength data will 
naturally overpredict the MIPS 24/tm flux, since we have 
forced a Type II template to fit the IRAC data. At the 
same time, this difference is not significant for most of the 
sample. The range of SPS+AGN templates that fit the 
observed UG7ZJK+IRAC data within the photometric 
errors is within ler of the MIPS uncertainty for the major- 
ity (five out of the eight) of our objects. Since we do not 
use the MIPS data in our SED fitting, the discrepancies 
between predicted and observed fluxes suggest that the 
Lboi values reported in Table [2] may slightly overpredict 
the true values. 

4.2. Correcting for AGN Emission 

To understand the true AGN host galaxy stellar popu- 
lations, it is very important to correct for AGN emission 
in the SED of those objects without IRAC data, where 
we can't directly quantify the AGN contribution. Many 
authors fit AGN host galaxy SEDs using SPS models 
alone (see iBundv et all 12008 iBluck et al.ll20H for dis- 
cussions), and find that while the presence of an AGN 
does affect stellar mass estimates, it is not a significant 
difference. Here we adopt a different approach, and use 
our results from the SPS+AGN modeling to quantify 
how much the presence of an AGN affects the stellar 
modeling parameters. For the 11 objects in the IRAC 
AGN sample, we refit the SED but only through the 
if -band, and using an SPS-only model. The difference 
in best-fit parameters based on if -band and SPS-only 
vs. if-band+IRAC and SPS+AGN modeling indicates 
how the majority of the objects in our sample without 
IRAC coverage must be corrected. It is valid to apply the 
correction factors derived from the IRAC AGN sample 
to our entire sample, given that the objects with IRAC 
coverage have similar average 1Z magnitudes and 1Z—K 
colors to those without IRAC coverage. The procedure, 



star-formation models, and parameters are the same as 
those discussed in ^3.21 but here we do not include the 
AGN template in the if -band-only fits. 

The resulting stellar masses and SFRs calculated with 
an AGN component are lower than those calculated with- 
out the AGN component. The best-fit ages and values 
for E(B — V)sf are consistent between fits with and 
without the AGN component. For a few objects, such as 
Q1623-BX663, the best-fit parameters are quite different, 
as this galaxy suffers from the most AGN contamina- 
tion at the shorter wavelengths (Fig. [6]). The results for 
HDF-BMZ1156 and HDF-BMZ1384 are also discrepant, 
due to the bright K magnitude for these objects, which 
is not well-fit in the dual-component modeling extending 
to IRAC wavelengths. 

We analyze the AGN contamination by calculating the 
ratio between the parameters (A/*, SFR, E(B — V)sf, 
and age) resulting from fits with the SPS+AGN modeling 
to the average parameters resulting from the SPS-only 
modeling, which we refer to as "correction factors." We 
calculated the ratios for the best-fit parameters, and wc 
also analyze the AGN contamination using results from 
the Monte Carlo error analysis ( §3.4j) . In this method, we 
calculated the average value for the full distribution of 
Monte Carlo output parameters, and then divided these 
averages for our SPS+AGN modeling by the averages 
from the SPS-only modeling. This second method gives 
a better understanding of the full set of models that fit 
an object's SED within the error on the photometry and 
therefore we only report correction factors calculated in 
this manner. We note however that the correction factors 
do not depend significantly on the method used. The 
correction factors calculated for the 11 IRAC AGNs are 
shown in Table |U While our results suggest that the 
UG1ZJK photometry for Type II AGNs is largely free 
from AGN emission, the best-fit masses and SFR values 
derived without including an AGN component are larger 
by a factor of 1.4 (although, on an object-by-object basis, 
these factors for mass and SFR are not always the same) . 
The values for E(B — V)sf are only larger by a factor of 
1.1, and the ages are consistent. We use these factors to 
correct the best-fit parameters based on SED modeling 
results for those objects where we do not have IRAC data 
under the assumption that the SPS+AGN model is more 
appropriate. 

We also calculated AGN correction factors using the 
CB12 and Mar05 modeling. The correction values calcu- 
lated for both the CB12 and Mar05 modeling are similar 
to those calculated using the BC03 modeling. While the 
mass and age correction values are smaller than than 
those derived from the BC03 models, they differ by less 
than one standard deviation of the distribution of BC03 
correction factors. The corrected best-fit parameters are 
not significantly dependent on the choice of model (e.g., 
BC03, CB12, or Mar05). 

4.3. Full Sample Modeling Results 

To fit the objects in our AGN sample lacking multi- 
ple IRAC photometric measurements, we used SPS-only 
models (using the same procedure described in §3.2[ but 
without the addition of an AGN model). For each of 
these objects, we corrected the best-fit mass, SFR, age, 
and E(B — V)sf for the presence of an AGN using the 
average correction factors reported in the previous sec- 
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Fig. 7. — SED models for the 17 objects that do not have IRAC data, or have IRAC data that is not sufficient for the purpose of discerning 
the presence of a power-law in the rest-frame near-IR. The models shown here are BC03 models, and do no t inc lude any contribution from 
an AGN. The best-fit mass for each object is given, corrected for the presence of an AGN as described in i|4,2l 



tion, and list the best fit values (both uncorrected and 
corrected) in Table [5] The errors on the corrected pa- 
rameters reflect both the errors derived using the Monte 
Carlo bootstrapping described in t )3.4l as well as the stan- 
dard deviation of the correction factors, given in Table 
[4j Objects with only three photometric data points have 
reduced x 2 values given by a dash. SED fits for each of 
these objects are shown in Figure [7] The masses indi- 
cated in these figures are the corrected best-fit masses. 
The x 2 values indicate that many of the objects are well 
fit by the SPS models. In the case of Q0000-C7 and 
Q1623-BX747, the SEDs were not well fit by the CSF 
models. For Q0000-C7, a blue G—1Z color leads to mod- 
els that are unable simultaneously to fit both the UV and 
near-IR portions of the SED. For Q1623-BX747, the J, 



K, and Channel 1 photometric datg0 show a decline to 
larger wavelengths, which is not reflected in the best-fit 
stellar models. 

We summarize the full statistics for our sample of UV- 
selected AGNs in Figure iJ All 28 of the AGNs with 
fitted SEDs are represented in these histograms. The 
AGNs with IRAC coverage are plotted with their best- 
fit parameters from the dual-component modeling, and 
the AGNs without IRAC coverage are plotted with the 
corrected best-fit parameters, calculated as described in 
£14.31 The average values for the best-fit parameters and 
standard deviations are indicated on the figure with ver- 



8 The data set for this object has coverage in one IRAC C hann el, 
which is not enough to determine an AGN power-law; see £|3. 1 1 
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Fig. 8. — Histograms of best-fit stellar parameters for our sample of UV-Selected AGNs. The histogram represents the total number 
of UV-selected AGNs, and we highlight the IRAC AGNs in red. Vertical solid lines indicate the averages for each parameter, and the 
dashed lines indicate the standard deviations of the distributions. The average (median) values for the best-fit parameters from all 
of the AGNs are, <log(M„/M Q )} = 10.85(10.71) ± 0.36, (SFR) = 63(37) ± 67 M Q yr" 1 , (log(Age/Myr)) = 3.19(3.18) ± 0.26, and 
(E(B - V)sf> = 0.22 (0.23) ± 0.12. 

tical lines. There is a large spread in E(B — V)sf 
values, with an average (median) of (E(B — V)sf) = 
0.22 (0.23) ±0.11, and a range of 0.00-0.40. The best-fit 
star-formation rates, with a mean of (SFR) = 63 (37) ±67 
Mq yr -1 similarly span a large range of values between 
5.2 - 309 M Q yr -1 . The most striking results, however, 
are the ages and masses for this sample. As suggested by 
the 71— K colors shown in Figure the best-fit masses 
are high: log((M*/M Q )) = 10.85 (10.71) ± 0.36. The 



AGNs that comprise the sample have old ages, with an 
average age of log((age/Myr}) = 3.19 (3.18)±0.26. These 
average stellar masses and ages are large compared to the 
typi cal values quoted for UV-selected non -active galax- 
ies (jErb et all [200611 iShaplev et al|[200l . and will be 
discussed further in $5] and Sjll 

We can use the L\, i values derived from the dual- 
component fitting to estimate accretion rates for the 
IRAC AGNs. It should be noted that, as described in 
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Fig. 9. — Observed AGN fraction as a function of stellar mass. 
The fractions were calculated by comparing the number of AGNs 
in a bin to the number of non-active galaxies in the full comparison 
sample. The AGN fraction rises steadily from log(M*/M0) = 10 
to log (M*/Mq) = 11.5, where it reaches a maximum of 19%. 
Uncertainties in observed AGN fraction were calculated by varying 
the masses of the AGNs by their mass errors and recalculating the 
AGN fractions for each bin. Horizontal error bars indicate the 
width of the bins. As there are AGNs out of 4 objects in the 
highest-mass bin, we estimate the 68% confidence limit on the AGN 
fraction (0.25) using binomial statistics. At the 95% confidence 
level, this fraction is 0.53. 

Section 14. 11 the values of L& ; for our sample are po- 
tentially too high as a result of the overprediction of 
the observed MIPS flux, so the derived accretion rates 
and Eddington ratios should be viewed as upper lim- 
its. The accretion rate is described by the equation 
M = Lb i/e ra dC 2 , where we assume a standard accre- 
tion efficiency e ra d = 0.1. Using the values for Lboi from 
our fitting, we calculate an average (median) accretion 
rate for our sample of 0.3(0.1) M Q yr _1 with a range 
of 0.04 — 0.9 M yr -1 . These rates are similar to the 
range of accretion rates fou nd for optically select ed Type 
I AGNs at 1 < z < 2.2 from lMerloni et all (j2010T ) (with a 
median accretion rate of M ~ 0.4M Q yr _1 ), and a factor 
of several lower th an those estimated fo r X-ray selected 
Type II QSOs from lMainieri et al.l (|201lD (with a median 
accretion rate of M ~ 1 M Q yr -1 ). 

The results from the SED fitting were also used to 
derive the properties of the central active galactic nu- 
cleus in each of our galaxies. Black hole masses are typ- 
ically estimated in Type I AGNs using broad emission 
line widths and AGN continuum luminosities by setting 
Mbh oc v 2 r, where r, the size of the AGN broad line re- 
gion, is estimated from the continuum luminosity using 
the relationship between luminosity a nd BLR size based 
on reverberation m apping studies (e.g. JBentz et al.f2009t 
iBarth et al.|[20ll . For our sample of Type II AGNs, 
black hole masses we re estimated using the Mbh — -^buige 
scaling relation from lHaring fc Rix (200J): 



\og(M BH /M e ) = (8.20±0.10)+(1.12±0.06)log(M 6tli3e /10 1 

(2) 

For this calculation, we assumed that the stellar masses 
derived from our SED fitting roughly corresponded to 
the bulge masses for these objects. We adjusted the 
black hole masses for redshift evolution following the 
param eterization der i ved fo r Type I AGNs out to z = 
2.2 in iMerloni et al.l (120101) (but see lLauer et al.l l2007t 
I Alexander et al.H2008alh 

Alog(M BH /M*)(z) = (0.68 ± 0.12) log(l + z) (3) 

The values for \og(M B H) for our sample of AGNs 
have a range of 7.7 — 8.7, with an average (median) 
of log(M Bff /M Q )= 8.36(8.29). The typical uncertainty 
in M BH is ~ 36%. For the eleven AGNs with IRAC 
data and Lhoi measurements, we calculated Eddington 
ratios (A E dd = Lboi/LEdd, where L E dd/£© = 3.2 x 
10 4 Mbh/Mq j3- The range is low er than that presented 
for X-ray selected Type II QSOs in lMainieri et al.l (|2011[ ) 
(in which A E dd values were calculated in the same man- 
ner). The median Eddington ratio for our sample is 0.03 
(with a typical uncertainty o f ~ 45%), whi le the me- 
dian Eddington ratio for the iMainieri et al.l sample is 
0.1. Therefore the black holes hosted by the AGNs in 
our sample are accreting at significantly sub-Eddington 
rates. 

5. AGN HOST GALAXY TRENDS 

In order to understand the impact of an AGN on its 
host galaxy, we must compare the UV-selected AGN host 
galaxy properties to those of a representative sample of 
non-active galaxies. As the AGNs studied here appear 
to be hosted by galaxies drawn from the same parent 
popu l ation as the non- AGN LBG sample (|Steidel et al.l 
l200l lAdelbereer et al.ll2005h . a comparison of the AGN 
and non-AGN host galaxy stellar populations in prin- 
ciple provides a means of testing the impact of AGN 
feedback. In this section, we examine a non-active sam- 
ple of UV-selected galaxies, within which we isolate a 
sample matched in stellar mass for a controlled compar- 
ison between AGNs and non- AGNs (i j5.1|) . Using these 
samples, we examine the relationship between SFR and 
stellar mass f £|5.2p . as well as rest-frame U — V color and 
stellar mass f ^5.3[) . Finally, making use of our unique 
spectroscopic data set, we search for connections between 
best-fit stellar population parameters and rest-frame UV 
spectroscopic features ( ij5.4[) . 

5.1. Non-Active Comparison Sample 

A careful examination of the demographics of the UV- 
selected AGNs requires the selection of a comparison 
sample of non-active galaxies. In Ej2j we highlighted a 
sample of non- AGNs that lie in the same redshift range. 
Recent results have demonstrated the importance of con- 
trolling for stellar mass in understanding the relationship 
betw e en AGNs and their host galaxies (jSilverman et al.l 
2009; IXue et al.ll2oTol: lAird et all 120121) . and so we also 

9 We calculated an Eddington ratio for Q1623-BX663 of 0.49, 
but this value is not meaningful due to the large calculated error 
on the stellar mass for this object. 



16 



HAINLINE ET AL. 



constructed a comparison sample that spans a similar 
mass range to that of the AGN host galaxies. For each 
AGN in our sample, we chose six non-active sources at 
random, without any duplicate objects, and character- 
ized by masses within the uncertainty of the selected 
AGN host mass. The size of the mass-matched sam- 
ple was limited by the nature of the stellar mass distri- 
bution of the parent non-AGN sample, in which there 
are only six unique objects for each AGN at the mas- 
sive end. Thus, the final mass-matched sample has 168 
objects, as compared to the sample of 28 AGNs with 
stellar population fits, and has similar redshift proper- 
ties to the AGN and non-AGN samples. The redshift 
distribution of the mass-matched sample is character- 
ized by (z) = 2.52 ± 0.37, which closely follows that 
of the UV-selected AGN sample. We will refer to the 
larger sample of non-AGNs as the "full non-active com- 
parison sample," and the subsample matched in stellar 
mass as the "mass-matched comparison sample." The 
average mass for the mass-matched comparison sample 
is log((M*)) = 10.85 ± 0.36, which, by construction, is 
the same as for the AGNs. We modeled the photome- 
try of non-AGNs using BC03 CSF SPS-only models as 
described in §4.31 

We used the full non-active comparison sample to cal- 
culate the observed fraction of AGNs as a function of 
stellar mass. In Figure |9l we show the results, with er- 
rors derived from the mass errors for the AGNs that went 
into each bin. These results indicate that the fraction 
of AGNs detected rises steadily as a function of stellar 
mass, peaking at 19% for log(M*/Af Q ) = 11.25. Above 
this mass, the number of non-AGNs in the full non-active 
sample drops to only 4 objects, preventing the determi- 
nation of tight constraints on the AGN fraction. We 
observe AGNs predominantly in higher mass galaxies, 
although as described in Sj6j this trend likely reflects our 
increasing incompleteness towards lower stellar masses. 
Accordingly, these results may indicate that at least 19% 
of UV-selected star-forming galaxies to 1Z— 25.5 host ac- 
tively accreting black holes. We return to this discussion 
in §B 

5.2. SFR - Stellar Mass Relation 

Star-forming galaxies follow a tight positive correlation 
between SFR and stellar m ass, which has been observed 
locally (IBrinchmann et"a . 



termediate (jNoeske et al 



2004 ), as well as out to in- 
2007TlElbaz et al l I2007T z 



1) an d high redshifts (iReddv et al l l2006at IDaddi et al.l 
[2007t iReddv et all 120121 z ~ ~ 3). This trend is 
used to argue for the smooth star-formation histories 
of the majority of s t ar-forming g alaxies in the universe 
(jNoeske et al.l l2007t IDaddi et al.l 120071 ). This relation 
is characterized by roughly the same sl ope across a 
wide r edshift range (although results from iNoeske et all 
(|2007l ) indicate a non-linear relation for star-forming 
galaxies at intermediate redshifts), while the normaliza- 
tion has decreased by an order of magnitude from z ~ 2, 
reflecting the overall dec rease in the global star-formation 
rate density with time (jDaddi et al.l 12007) . We plot the 
SFR and stellar mass for our AGNs with red points in 
Figure [TUl We find a positive correlation for the AGN 
host galaxies with a best-fit slope of 0.20±0.04, shallower 
than derived for non-AGNs. 
We compared the AGN population to our non-active 
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Fig. 10. — SFR, vs. stellar mass for the AGNs (red points) along 
with the full comparison sample of non-active star-forming U G1Z- 
selected galaxies (grey points) and the mass-matched comparison 
sample (blue points). AGNs are found at the high mass end of the 
full relation. However, both AGNs and mass-matched non-AGNs 
span a similar range of SFR values at a given stellar mass. The 
ridge of grey points at high SFRsis a result of imposing a lower 
age limit of t$F > 50 Myr (see M5 2JI . 



comparison samples in order to examine the origin of the 
shallow slope. In Figure [10] we plot log(SFR) against 
log(Af„.) for the non-active comparison sample in grey 
and for the mass-matched comparison sample in blue. 
The non-active galaxies display a positive trend with a 
large scatter. The ridge of points seen at the high SFR 
side of the relation is an artifact of imposing a minimum 
age of tsF = 50 Myr on our fits. Given the linear re- 
lation between M* and SFR for CSF models, all points 
with best-fit ages of t$F — 50 Myr will fall on a line 
of constant iW*/SFR. We note that the sharp cutoff for 
galaxies with low SFRs at a fixed stellar mass is due to 
the restriction on the maximum allowed age, which must 
be younger than the age of the universe at the redshift 
of each object. The lack of points at the low-mass, low- 
SFR end is due to a bias in the U GTZ selection technique 
towards selecting those galaxies with the highest SFR at 
a given mass, causing the slope of the relation to be come 
shallower at the low-mass end (|Reddv et al.|[2012T) . At 
the high-mass end, the shallow slope reported for the 
AGNs is similar to what is observed among the non- 
active sources. Above log(M*/M Q ) = 10.5, the UGTZ se- 
lection is biased against selecting those galaxies with the 
largest SFR values at a given stellar mass, as such objects 
will have large amounts of dust extinction and red U GTZ 
colors, placing them outside the U GTZ selection window. 
This bias tends to flatten the log(M») — log(SFR) rela- 
tion at the high-mass end. Of particular interest, the 
SFR values spanned by the UV-selected AGN sample 
(((SFR)) = 63 (37) ±67 M Q yr" 1 ) are similar to those in 
the mass-matched comparison sample, where the average 
(median) SFR is (SFR) = 45 (35) ± 38 M yr" 1 . These 
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Fig. 11. — Rest-frame U — V color vs. stellar mass. The rest- 
frame U — V c olors are calculated from the best-fit SEDs as de- 
scribed in i|5,3l AGNs are indicated in red, the non-active compar- 
ison sample in grey, and the mass-matched sample in blue. The 
objects plotted were initially selected by their observed U GR colors 
to be strongly star forming, which leads to the tight blue sequence 
IILabbe et al.|[2007T) . The UV-selected AGN hosts are located in a 
similar color-mass space to the mass-matched sample, but are ab- 
sent below log(Af*/M0)~ 10. In the inset, we plot the histograms 
of U — Vrest for the three samples, using the same colors. These 
histograms have been normalized to have the same peak values for 
clarity. 

similarities are observed over the full AGN and mass- 
matched comparison samples, as well as in bins of stellar 
mass. We quantify them using the Kolmogorov-Smirnov 
two-sample test, which yields a 44% probability that the 
AGN and mass-matched non-AGN SFR distributions are 
drawn from the same parent population. 

The above trends can also be considered in terms of the 
specific star-formation rate (sSFR = SFR/M*, which, for 
CSF models is equal to 1/tsF)- We calculated an average 
(median) sSFR = 1.36(0.66) Gyr" 1 for the AGN hosts, 
with a range of sSFR: 0.4 — 9.2 Gyr -1 , indicating active 
star-formatiorFI. We find that the sSFR values for our 
mass-matched non-active sample span a similar range, 
with an average (median) sSFR =1.10 (0.47) Gyr -1 , and 
a range of sSFR: 0.33 - 8.78 Gyr" 1 . The sSFR values 
for our Type II UV-selected AGNs are similar to those 
obtained for a sample of X-ray selected Type II QSOs 
with a ctive star-formation at z ~ 2 from M ainieri et al.l 
poTl (sSFR = 0.8-3 Gyr -1 for the high-mass objects at 
z ~ 2). These authors found consistent evolution for the 
sSFRs of the QSO host galaxies and thos e of non-active 
BzK-s elected star-forming galaxies from iPannella et al.l 
(2009). Similarly, once we control for stellar mass, the 
SFR and sSFRs of AGNs are indistinguishable from those 

10 It should be noted that the minimum and maximum allowed 
age limits for our SED fitting translate into limits on the range of 
sSFR values that can be recovered for both AGNs and non-AGNs. 
These lim its are apparent in the plot of SFR vs. M* shown in 
Figure [TO] 



of non-AGNs. 

5.3. Color - Stellar Mass Relation 

A basic probe of the demographics of AGN activity is 
provided by the optical color magnitude diagram, which 
can also be cast in terms of galaxy color as a function 
of stellar mass. In the local universe, galaxies display 
a bimodality in their observed colors, separating into a 
population of lower-mass blue star-forming galaxies and 
a sequence of high er-mass red quiescent galaxies (e.g., 
iBaldrv et all 120041 ). This same bimodality ha s been ob- 
serve d in galaxies at intermediate fa ~ 1. | Bell et al.l 
20041) and higher redshifts (z ~ 2. iCassata et al.l l20oT 



Kriek et al.l 120081: iWilliams et "all 12000) • The "existence 
of these two populations out to z ~ 2 suggests that 
the mechanisms driving the evolution of galaxies from 
blue to red began early in the history of the universe. 
AGNs have been proposed as a potential driver of the 
evolution in the color-mass relation by providing en- 
ergy to the ISM that quenches star-formation, causing 
a galaxy to tra nsition from the blue population of galax- 
ies to the red (jCroton et al.l 120061 : iHopkins et all 120081) . 
In support of this scenario, X-ray and optically selected 
AGNs at low and intermediate redshifts are found to re- 
side in galaxies in the transition region i n optical color 



space between the blue and red sequences dNandra et al 
200% I Salim et al.1 120071: ICoil et al.1 l200l IHickox^Tal 
200a ISchawinski et al.1 l2O10f T However, recent AGN 



host galaxy studies indicate that both mass selection 
and dust obscuration effects may bias c onclusions about 
the ro le of AGNs in galaxy evolution. iSilverman et al.l 
(2008) demonstrated that a sample of X-ray and mass 
selected AGNs were preferentially found in galaxies with 
ongoi ng star formation and blue optical colors. Simi- 
larly, iXue et al.l (poToh . using a mass-matched sample 
of X-ray selected AGNs, showed that the fraction of 
AGNs remains constant as a function of color. lAird et all 
(2012) found that, when stellar-mass selection effects are 
taken into account, AGN hosts are preferentially found 
in galaxies with blue and gree n optical colors, but can 
exist in galaxies of any color. iCardamone et all ()2010l ) 
corrected AGN host colors for the presence of dust, re- 
vealing that AGN hosts cleanly separate into a popula- 
tion of young, dusty, star-forming galaxies, and older, 
red, passive galaxies. 

We used the results from AGN SED modeling (gj) to 
examine the positions of the UV-selected AGNs in color- 
mass space. We calculated rest-frame colors of AGN host 
galaxies using the best-fit SPS models, along with the 
Johns on UBV transmission curves from lMaiz Apellanid 
(2006). The estimation of the rest- frame colors re- 
quires an assumption that the models accurately rep- 
resent the stellar populations in the rest-frame UV and 
optical, which is supported by the good SPS+AGN dual- 
component modeling fits to the IRAC AGNs. In order to 
compare the AGN host colors to those of the non-active 
UV-selected star-forming galaxies, we also calculated the 
rest-frame colors for the non-active comparison samples 
described in the previous section. The (U — V) res t,AB 
vs. M* relation is shown in Figure [HI where the AGNs 
are plotted in red and the non-active galaxies in grey. 
The non-AGN hosts follow a trend where the most mas- 
sive galaxies are also the reddest, similar to the blue se- 
quence observed for galaxies out to z ~ 3 (Labb l~et al.l 
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Fig. 12. — Histograms of best-fit stellar population parameters. 
Distributions for AGNs are shown in red, those for the full non- 
AGN sample in grey, and for the mass-matched sample in blue. 
These histograms have been normalized to have the same peak 
values for clarity. The AGNs have similar properties to those of 
the mass-matched sample on average. The AGNs have a slightly 
higher average extinction ({E(B — V),$f.AGn) = 0.22 ±0.11) than 
the mass-matched sample ({E(B — V)sf,mm) = 0.20 ± 0.09). 



2007). The sample of UV-sclectcd AGNs are preferen- 
tially hosted by high-mass galaxies and these objects are 
found to span the same range of (U — V) rest ,AB colors 
as those of the non-AGNs in the mass-matched sam- 
ple. The average (median) color for the AGN sample 
is (U - V) rest ,AB = 0.9 (0.9) ± 0.2, while for the mass- 
matched sample, it is (U_ — V) r e „ t.AR = 1-0 (1-0) ± 0-3 
(inset of Fi gure ITT1). Kue et all (|2010D . iMainieri et all 



(201 3), and lAird let all (|2012H found similar results for 
the demographics of their AGN host galaxy samples. 
All of these results demonstrate the importance of mass- 
selection in the comparison of AGN hosts to non-active 
galaxies. The redder colors of the AGN hosts relative 
to the full non-AGN comparison sample do not indi- 
cate a lack of ongoing star formation, but, rather, as we 
show below, more dust and perhaps older stellar popula- 
tions characteristic of the massive end of the UV-selected 
star-forming population. The UV-selected AGNs were 
initially chosen based on their UGH colors to be star- 
forming, and, as discussed in the previous section, we 
found that all of the AGNs have sSFR > 0.4 Gyr -1 , in- 
dicating active star formation. Figure [12] demonstrates 
that the SFR, age, color, and dust extinction proper- 
ties are similar between the AGN sample and the mass- 
matched comparison sample. The presence of an AGN 
does not seem to affect the global properties of their host 
galaxies in comparison to the mass-matched sample. The 
lack of UV-selected AGNs in host galaxies below 10 10 Af Q 
will be discussed in 

We can use the mass-matched sample to explore the 
7Z— K distribution shown in Figure[2]in more detail. The 
fact that AGNs occupy the red portion of this distribu- 



tion can be explained by their massive host galaxies, and 
the relation between rest-UV/optical color and stellar 
mass. The average 1Z—K color for the UV-selected AGN 
sample is 4.0 ± 0.7, while for the mass-matched sample, 
the average is 3.6 ± 0.6. We used the results from our 
dual-component modeling to estimate that the presence 
of an AGN caused the 1Z—K color to be 0.3 mag redder 
than it would be in the absence of an AGN. If we sub- 
tract this difference from those AGNs where we did not 
use dual-component modeling, the average UV-selected 
AGN host galaxy 7Z—K color becomes 3.6 ± 0.7. There- 
fore, at fixed stellar mass, the AGN host galaxies are 
typical in terms of their UV/optical (i.e., 1Z—K) colors. 

5.4. UV Composite Spectra 

In lHainline et all (|2011a[) . we investigated the UV 
spectral properties for our sample of UV-selected AGNs. 
The UV portion of the galaxy spectrum contains strong 
emission and absorption features that we used to probe 
both AGN and outflow activity. At z ~ 2 — 3, the indi- 
vidual UV spectra have low S/N, and, to overcome this 
limitation, we averaged the individual spectra to create 
higher S/N composite spectra. One of the most striking 
results from lHainline et all (|2011al) consisted of the ex- 
tremely red rest-frame UV continuum of the AGN com- 
posite spectrum, relative to that of the z ~ 3 non-AGN 
composite spectrum from iShaplev et all (|2003f ). We can 
now use the results from SEE) modeling to investigate 
the physical origin of this difference. Furthermore, we 
can analyze how the rest-frame UV spectroscopic prop- 
erties of AGN vary as a function of stellar population 
properties. 

For this analysis, we created composite spectra from 
sub-samples of AGNs separated by the host galaxy prop- 
erties. Due to the small sample size of objects with SED 
fitting parameters (28 objects) and our desire to maxi- 
mize the S/N of the resulting subsample composite spec- 
tra, we simply divided the sample in half for these anal- 
yses. We const ructed composi t e spec tra following the 
methodology of lHainline et all (|2011al) . Each individ- 
ual spectrum was first shifted to the rest frame using an 
estimate of the systemic redshift from He IIA1640. We 
converted each spectrum from the units of flux density to 
those of luminosity density, scaled each one to a common 
median in the wavelength range of 1250 - 1380 A, and 
combined the spectra, with the four highest and lowest 
outliers removed at each wavelength. 

In order to investigate the red UV continuum of the 
AGN composite spectrum, we consider the properties 
of the mass-matched comparison sample. Using exist- 
ing rest-frame UV spectroscopy for the full non-active 
sample, we created a UV composite spectrum from the 
objects in the mass-matched comparison sample. We 
used the same method to create this non-AGN com- 
posite as for the AGN composite. In Figure [T31 we 
compare the full AGN composite spectrum to the mass- 
matched non-AGN composite as well as the non-AGN 
LBG composite at z ~ 3 from IShaplev et all ((2003). 
The UV continuum shape is commonly described by /?, 
the slope of a power law of the form L\ oc A^, which 
is fit to the continuu m. For the AGN com posite spec- 
trum, (3 = -0.3 ± 0.2 (|Hainline et alll2011al) . which is a 
much redder s l ope co mpared to the 1BG composite from 
IShaplev et all (l200l . where (3 = -1.5. We measure a 
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Fig. 13. — The UV composite spectra for the UV-selected active a nd non-active samples . We plot the AGN composite in red, the mass- 
matched composite in blue, and the non-AGN LBG composite from Shaplcv et al. (2003) in grey. The power-law slope, ft, as measured 
directly from the AGN spectrum, is /3 = —0.3 ±0.2, while for the mass-matched composite, the slope is = —0.7 ±0.1. The similar slopes 
may be due to the similar levels o f dust extinction in t hese two samples. The mass-matched composite shows a much redder slope than 
the full non-active composite from Shaplcv ct al. (2003), which has a measured power-law slope of /3 = —1.49. These spectra confirm our 
conclusion from the photometry that the AGN UV continuum is primarily due to stellar light, rather than AGN emission. 

value of /3 = —0.7 ±0.1 for the mass-matched composite, 
redder than the full non-AGN compositfl. The com- 
parably red slopes of the AGN and mass-matched UV 
composite spectra (in contrast to the sig nificantly bluer 
slope of the total LBG composite from iShaplev et al.l 
2003) confirms that the underlying UV continuum seen 
for the AGN sample is predominantly due to stellar 
light, and not AGN emission. Furthermore, since scat- 
tered light from a buried AGN would have an intrinsi- 
cally bluer rest- frame UV spectrum t han that of a stel- 
lar population (Zakam ska et al.ll2006D . the slight differ- 
ence in UV slopes can not be attributed to AGN emis- 
sion. We can naturally explain the redder slopes ob- 
served in the high-mass active and inactive samples by 
looking at the distributions of reddening for each sub- 
sample. In Figure (T2J we show histograms of the best-fit 
properties of the AGN sample along with the full non- 
AGN and mass-matched comparison samples. As dis- 
cussed in g31 (E(B - V) S f,agn) = 0.22 ± 0.11 for the 
AGNs, while for the mass-matched comparison sample, 
(E(B - V) S f,mm) = 0.20 ±0.09. For the sample of non- 
AGNs at z ~ 3 from IShaplev et al.1 (|2003l ). the estimated 
(E{B - V) S F,Shapiey) = 0.15 ± 0.00 The errors on 
these values represent the standard deviation of the dis- 
tribution of values. The average extinction values for the 
AGN and mass-matched samples are signi ficantly differ- 
ent from the value for the non-AGNs in ShaDlev et al. 



1 Errors on the f} values were calculated following a bootstrap 
technique where 500 fake composite spectra were constructed from 
the sample of actual spectra and we measured ft for each fake com- 
posite. The error represents the standard deviation of the distri- 
bution of measured /3 values 

12 These E(B — V)sf value s were calculated dire ctly from the 
UV photometry as described in SEaple y et al. (2003). 



( 2003) , reflecting what is observed for the UV power- law 
slopes. 

We further investigated the role that extinction plays 
in the slope of the UV spectrum by separating our ob- 
jects into two bins based on E(B — V)sf- The AGNs 
were separated at the median extinction from our fits, 
E(B — V)sf = 0.242, and both composites are plot- 
ted in Figure HU In red, we plot the composite spec- 
trum for those objects with E(B — V)sf > 0.242, and 
in black, the composite spectrum for those objects with 
E(B — V)sf < 0.242. The composite spectrum for ob- 
jects with higher E(B — V)sf is much redder, with (3 = 
0.0±0.3. For the AGNs with lower E(B-V) SF , the com- 
posite spectrum has a (3 = —1.0 ± 0.3. This trend is also 
seen in the mass-matched comparison sample separated 
by E(B — V)sf, where the high-extinction mass-matched 
composite is redder, with /? = —0.4 ± 0.2, while the low- 
extinction mass-matched composite has /? = — 1.1 ± 0.1. 
The slope of the low-extinction mass-matched compos- 
ite spectrum {{E(B - V)sf,mm) = 0.12 ± 0.05) is sim- 
ilar to the slope of the low-extinction AGN composite 
{{E(B -V) SF ,agn) = 0.13 ±0.08), which provides more 
evidence that dust extinction significantly modulates the 
UV power-law slope in our sample of AGNs. 

We also consider how rest-frame UV AGN spectro- 
scopic properties vary with stellar mass. The continuum 
normalized composite spectra for t he A GNs separated 
by stellar mass are shown in Figure II 4 We measured 
the luminosities and EW for the strongest UV emission 
lines (H I Lya, N VA1240, C IVA1549, and He IIA1640) 

3 Given the correlation between stellar mass and E(B — V)sf, 
a separation by stellar mass divides the sample in almost the same 
way as a separation by E(B — V)$f- 
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Fig. 14. — Comparison of the continuum slope between the AGN composite spectrum for objects separated by E(B — V)sf- The composite 
with larger E(B — V)gp is significantly redder than the smaller E(B — V)sf composite. The UV power-law slope /3, as measured directly 
from the E(B - V) S F > 0.242 spectrum, is = 0.0 ± 0.3, while for the E(B - V) S F < 0.242 spectrum, /3 = -1.0 ± 0.3. The large average 
E(B — V)sf for our AGN sample indicated in Figure [8] leads to the red UV continuum seen in AGN UV composite spectrum. 



in both the high-mass and low-mass composites. Uncer- 
tainties on these values were calculated following a boot- 
strap technique where 500 fake composite spectra were 
constructed from the sample of sp ectra used in creat- 
ing the real composit e spectra (see iShaplev et all 120031 : 
lHainline et al.ll2011a|) . The results are shown in Table 
El We find that the EW values for the C IV and He II 
emission lines are stronger in the high-mass composite 
than the low mass composite, while the N V and Lya 
lines have statistically equivalent EW values. The EW 
of an emission line was calculated by integrating the lu- 
minosity in the line and dividing by the average of the 
luminosity density of the continuum on either side of the 
line. As these are Type II AGNs, the observed UV con- 
tinuum arises from stellar emission, which is supported 
by the similar UV slopes seen in the AGN UV composite 
and the mass-matched UV composite. We find that the 
emission lines used to identify the presence of an AGN 
(C IV and He II) are stronger in the high-mass compos- 
ite spectrum. This trend of decreasing high-ionization 
EW with decreasing stellar mass may explain the dearth 
of low-mass UV-selected AGNs in our sample. We will 
explore the origins and implications for this trend in the 
following section. 

6. THE DEMOGRAPHICS OF UV-SELECTED AGNS AT 
Z ~ 2 - 3 

We find an apparent segregation of AGN hosts at the 
high-mass end of the z ~ 2 — 3 galaxy population, as 
shown in Figures [TU1 [TTJ and [T3J At lower-redshift 
(z ~ 0.2 — 1.0), the preferential incidence of X-ray- 
sel ected AGN s in hi gher- mass host galaxies is explained 
bv lAird et all (pOlU ) as a selection effect, arising from a 
constant Eddington-ratio distribution and an X-ray flux 
limit. For a fixed Eddington ratio, AGNs in higher-mass 
galaxies have higher AGN luminosities (assuming a linear 
correlation between Mbh and M^ ost ), and are therefore 
easier to detect. We use the extinction-corrected C IV 



luminosity as a proxy for AGN luminosity. Without di- 
rect X-ray or UV-optical continuum measurements, we 
must resort to less direct measures of luminosity for our 
obscured sample. Here we assume that the extinction- 
corrected C IV luminosity from the narrow-line region 
traces the accretion rate onto the black hole, as shown 
for narrow AG N o ptical emission lines (e.g., [O III] A5007) 
by [Yel (|1980f ) and lShuderl (|1981| ). With this as s umpt ion, 
we can apply the arguments from lAird et all ()2012h to 
our higher-redshift AGN sample, in order to investigate 
the cause of the apparent segregation in host galaxy stel- 
lar mass, and the detectability of AGN activity in lower- 
mass host galaxies. The lack of detected AGN activity in 
the low-mass UV-selected sample may be explained sim- 
ply as a similar selection effect. That is, if black holes 
are smaller on average in low-mass galaxies, then their 
emission may fall below our detection threshold. 

In the first place, we must establish empirically how 
emission-line EW and luminosity vary with stellar mass. 
The EW of an emission feature depends on both the in- 
trinsic luminosity of the feature as well as the contin- 
uum luminosity. In Figure IT6l we plot emission line EW 
against stellar mass for objects in our AGN sample. The 
N V, C IV, and He II emission lines were measured from 
individual spectra. In the top panels, we also plot the 
EW values measured from the composite spectra sepa- 
rated by mass. For the luminosity values, the individ- 
ual spectra were flux-normalized such that the 1500 A 
continuum luminosity matched the estimated rest-frame 
1500 A luminosity calculated from G and/or 7£-band 
photometry for these objects. We then corrected each 
spectrum for extinction under the assumption that the 
stellar extinction estimated from the SED modeling ap- 
plies to the narrow-line region. The N V and C IV emis- 
sion line luminosities shown in Figure [16] demonstrate a 
weak trend with mass (Spearman's rank coefficients of 
0.45 and 0.43 for the N V and C IV distributions re- 
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Fig. 15. — Comparison of the strong emission lines between the AGN composite spectrum for objects separated by M*. These spectra are 
continuum normalized to better illustrate the line EW differences. The high-mass composite spectrum is shown in red, while the low-mass 
composite spectrum is shown in black. The AGN emission lines of N V, C IV, He II are stronger in the high-mass composite, with specific 
EW value given in Table [6] 

at \og(M if /MQ) ~ 11 to compute the average line lumi- 
nosity for two bins ol stellar mass below the minimum 
mass in our AGN sample: log(M*/M ) = 9.0 — 9.5 and 
log(M*/M Q ) =9.5-10.0. Based on the results of the 
SED modeling for the objects in the non-AGN compar- 
ison sample in these two stellar mass bins, we applied 
an average level of extinction ((E(B — V)sf) = 0.10 for 
the lower mass bin, and (E(B — V)sf) = 0.12 for the 
higher mass bin) to predict extinction corrected C IV 
luminosities. Finally, we used the objects in the lower- 
mass bins from the non-AGN comparison sample to esti- 
mate an average rest-frame 1500 A continuum luminosity 
density from the G and 1Z photometry for the objects. 
Based on the predicted C IV luminosities and typical 
UV continuum luminosity densities, we estimated the 
typical C IV EWs for objects in these lower-mass bins. 
The resulting EW values are EWciv = 3.4 A for the 
log(M*/M Q ) = 9.0 - 9.5 bin, and EW C iv = 6.6 A for 
the log(M*/M Q ) = 9.5 - 10.0 bin. 

With average C IV EWs in hand, we can now ask 
whether we would have identified these lines in galax- 
ies with log(M*/M Q ) < 10. To test whether these lines 
would be detectable in lower-mass galaxies, we used the 
individual UV spectra for 15 representative non-AGNs 
in each of these lower-mass bins. We calculated the un- 
certainty on the EW in the spectral regions around the 
C IV line using the rms noise in the individual spectra, 
and estimated that in order for the line to be detected, 
EWciv > 8.5 A for the log(M*/M ) = 9.0-9.5 bin, and 
EWciv > 9.0 A for the log(M*/M ) = 9.5 - 10.0 bin. 
These estimates represent a 5cr detection limit that was 
used for identifying C IV in the AGN spectra. For typical 



spectively, corresponding to ~ 2a significance) such that 
the higher mass galaxies have higher luminosity emission 
lineS 

Of the three emission lines used to select an AGN, 
C IV is typically the strongest, and we focus on that line 
in order to test the detectabili ty of lower-mass A GNs. As 
discussed above, results from lAird et al.l (|2012[ ) suggest 
that X-ray selected AGNs out to z ~ 1 have a univer- 
sal distribution of Eddington ratios independent of host 
galaxy stellar mass. If AGNs at z ~ 2 — 3 are also charac- 
terized by a universal Eddington ratio distribution, and, 
furthermore, if Mbh and M* are correlated as in the 
low-redshift universe, galaxies with higher stellar masses 
should exhibit AGN-related emission lines with larger 
luminosities (i.e., a constant ratio of Lciv/M*). The as- 
sumption of a constant Lciv/M* ratio is supported by 
the observed correlation between C IV luminosity and 
stellar mass in our sample. Using this assumption along 
with the assumption of the correlation between Mbh and 
M„, we can test our ability to detect UV-selected AGNs 
in host galaxies with a given log(M*/M Q ). 

Since AGNs are easiest to find at the highest stellar 
masses, we use the spectra of the highest-mass AGNs in 
our sample to estimate the ratio of Lciv/M*. We av- 
eraged the C IV line luminosities for all AGNs in our 
sample with log(M»/M Q ) > 11. Under the above as- 
sumption that emission line luminosity scales linearly 
with stellar mass, we used the average line luminosity 

14 This trend is strengthened considerably by removing Q2233- 
MD21, which has very luminous emission lines and a low best-fit 
stellar mass. The Spearman rank coefficients rise to 0.68 and 0.62, 
for the N V and C IV distributions respectively, corresponding to 
~ 3cr significance. 
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Fig. 16. — EW (top three panels) and line luminosities (bottom three panels) for the strong UV emission lines in the AGN sample plotted 
against stellar mass. The black points are measured from the individual AGN spectra, and the red points in the EW plots represent values 
measured from the composite spectra as separated into two mass bins. For C IV and He II, the points from the composite spectra indicate 
that the EW values are higher in the high mass composites. The luminosity values are corrected for dust extinction using the E(B — V)sf 
values from the SED fitting, and show a positive trend with mass, such that the higher mass galaxies have intrinsically brighter emission 
lines. 

spectra of objects at log(M H ,/Af ) < 10.0, the predicted 
C IV EW values would therefore be too low to be de- 
tected. In conclusion, if galaxies with log(M*/M ) < 
10.0 hosted black holes accreting at the same average 
Eddington ratio as in higher mass galaxies, we would be 
unable to detect them given the emission line sensitivity 
in our spectroscopic sample. 

The segregation of UV-selected AGNs in high-mass 
(log M*/M Q > 10.0) host galaxies is naturally explained 
based on two key assumptions: (1) that {Lew I Mbh) 
is a constant as a function of Mbh (i-e., a constant Ed- 
dington ratio distributi on), as demonstr ated for X-ray- 
selected AGNs at z < 1 (|Aird et al. 2012) and supported 
by the observed correlation between Lciv and M* in our 
sample, and (2) that Mbh is proportional to M*. Even if 
lower-mass AGN host galaxies were characterized by sim- 
ilar Eddington ratios to those at higher-mass, we would 
be unable to detect them based on their high-ionization 
rest-frame UV emission lines. Viewed from a different 
perspective, the very segregation of UV-selectcd AGNs 
in high-mass hosts suggests that Mbh and M* are al- 



ready correlated at z > 2, during the epoch when both 
bulges and black holes are actively growing. Alterna- 
tively, our results can be explained if there are no black 
holes at z ~ 2 — 3 in galaxies at log(M*/M©) < 10, or 
if the average Eddington ratio decreases in less massive 
systems. This latter s cenario is inconsis tent with results 
at lower redshift from lAird et al.l (|2012f ). 

The discussion thus far has focused enti r ely o n 
Type II AGNs. As present ed in iSteidel et al.l (2 002). 
Adclb erger fc Steidell |2005an . and lAdelberger fc Steidel 
(2005b]), within the UV-selected galaxy population, 
broad-line or "Type I" AGNs have a similar frequency 
to that of the Type II AGNs discussed here. It is im- 
portant to consider if any of our conclusions regard- 
ing the demographics of AGN activity at z ~ 2 — 3 
are impacted by including a description of these Type 
I AGNs. Although it is not possible to make robust 
models of the host galaxy stellar populations of Type 
I objects, due to significant contamination from the ac- 
tive nucleus, we can investigate the typical Eddington 
ratios for Type I objects spanning the same range in 
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Lf, i as that of our Type II sample. We estimated L^oi 
for the Type I objects from their rest-frame 1350 A lumi- 
nosities (i.e., ALa(1350A), calculated using the observed 
G- and 7?.-band photometry for these objects), which we 
then mu ltiplied by the bolometric c orrection factor of 
4.3 from IVestcrgaa rd fc Osmeil (|2009Q . We used the re- 
lationship between Lb d and Mbh implied b y the cor- 
relations presented in [Adclbcrg er fc Steidell (|2005bl ) to 
estimate Mbh values. We find typical black hole masses 
of \og(M BH /M Q ) ~ 8 - comparable to the black hole 
masses presented section 4.3. Accordingly, the Type I 
objects in the a matched range of L^oi must have similar 
Eddington ratios to those of the Type II objects pre- 
sented here, and appear to represent their unobscured 
analogs. It is worth noting that the L^oi distribution of 
Type I AGNs identified within the the UV-selected sur- 
vey extends to higher luminosities than those spanned 
by the Type II objects, with roughly one third of Type I 
sample characterized by Lboi > 10 46 erg s _1 . Eddington- 
limited accretion in these systems corresponds to black 
hole masses of log(MsH /Mq) ~ 8.4, which represents a 
lower limit in black hole mass. Accordingly, it seems un- 
likely that these black holes are hosted by significantly 
lower-mass galaxies than the Type II AGNs discussed 
in this work. It is remarkable that both types of UV- 
selected AGNs at fixed bolometric luminosity are charac- 
terized by Eddington ratios < 0.1 for typical black hole 
masses of \og(MsH /Mq) ~ 8. At some earlier time, 
these black holes must have been accreting at signifi- 
cantly higher Eddington ratios. An important goal will 
be the identification of the more active progenitors of the 
UV-selected AGNs at z ~ 2 - 3. 

7. CONCLUSIONS 

Using SPS modeling, we have examined the host 
galaxy populations for a sample of z ~ 2 — 3 Type- 
II AGNs identified by the presence of rest-frame UV 
emission features. The results of our AGN+SPS dual- 
component modeling demonstrate the importance of us- 
ing long wavelength coverage to properly model the SEDs 
of type II AGNs. Without data longwards of 1 — 2{im 
(rest frame), the masses and SFRs for our sample of 
AGNs would be overpredicted by an average factor of 
1.4. We have used the results of our modeling to exam- 
ine the host galaxy trends for our sample. Our primary 
results are the following: 

1. The host galaxies for the UV-selected AGNs 
have very high average (median) masses 
(log((M*/M )) = 10.85(10.71) ± 0.36) and 
high SFRs ((SFR) = 63 (37) ± 67 M yr^ 1 ). 

2. We have constructed a mass-matched non-AGN 
sample with which to compare to the AGN hosts. 
The AGNs lie at the high-mass edge of the U — 
Vrest,AB vs. M* relation compared to the full non- 
AGN sample, but span the same range of color as 
the mass-matched non-AGN sample. Additionally, 
the AGNs have similar SFRs to the mass-matched 
sample, indicating that the presence of an AGN 
does not affect the colors or star-formation prop- 
erties of these highly star-forming hosts at these 



redshifts. 

3. A primary result from Mainline et al.l (|2011al) is 
that the UV composite spectrum for the AGNs has 
a redder c ontinuum than the f ull non-AGN com- 
posite from[ShiiIixiLal. (2003). TheUV compos- 
ite constructed for the objects in the mass-matched 
comparison sample has a red continuum that is 
very similar to that of the AGN composite, and 
we argue that this result is due to more dust ex- 
tinction in the higher-mass galaxies hosting AGNs 
as well as the mass-matched sample. This similar- 
ity between the UV continua of the mass-matched 
non-AGNs and AGNs confirms that the UV con- 
tinuum of AGN hosts is dominated by starlight. 

4. We separated our AGNs by stellar mass into two 
subsamples and created composite spectra for each 
of these subsamples. The C IV and He II EWs are 
larger in the high-mass composite spectrum than in 
the low-mass composite spectrum. This trend pro- 
vides evidence that the EWs of the UV emission 
lines used to identify the AGNs depends on stellar 
mass. We also observe that the dust-corrected lu- 
minosities for the N V, C IV, and He II emission 
lines scale with host galaxy stellar mass. These 
trends would be predicted if the mass of the host 
galaxy scaled with the black hole mass and emis- 
sion line luminosity. 

5. We estimate the C IV line EW for potential AGNs 
hosted by galaxies of mass log(M*/M Q ) = 9.0 — 9.5 
and log(M*/M Q ) = 9.5 — 10.0, under the assump- 
tion that line luminosity scales linearly with galaxy 
mass. The predicted C IV EWs (3.4 A for the low- 
est mass bin, and 6.6 A for the higher mass bin) 
would not be detectable at the 5a level required for 
AGN identification in our sample, indicating that 
we cannot identify UV-selected AGNs in galaxies 
with stellar masses log(M*/M ) < 10.0. Addi- 
tionally, we can use the segregation of UV-selected 
AGN hosts to argue that Mbh and M* are already 
correlated at z > 2. 

While we have demonstrated that AGNs hosted by 
lower-mass galaxies would not be detected in our sample 
based on the presence of UV emission lines, AGNs can be 
identified using X-ray luminosity and IR colors. In future 
work, we will investigate the demographics of AGNs se- 
lected using these additional methods and obtain a more 
complete picture of AGN activity at z > 2. 
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08 


n 00 1 n 1 c 

0.88 ± 0.16 





32 


_l_ 

zt 


u. 


07 
.0 1 




no.no. c: 1 c 
02:0o:5o.lo 


1 1 . op . on 1 
ll:oo:o0.1 


2.357 


25. 


n 1 


4- 

Zt 


0. 


34 







56 


-1- 
ztz 


n 

u. 




/~\ nor/"' tv r t~\ ^7 

Q0256-MD37 


no cr n no 01 

02:59:02.21 


nn 1 no a 

00:12:03.4 


2.803 


24 


oo 


_l_ 

Zt 





18 







63 


-1- 

ZL 


0. 


1 A 


/~\nnoo n /rnoo 

O0933-MD38 


nn 00 /i c\ 

09:33:48.00 


OO A A OO O 

28:44:32.3 


2.763 


22 


.01 


1 

Zt 


0. 


06 







04 


_l_ 

zt 


n 
u. 


no 
.uy 


(41217-BX4o 


TO T n 1 t\ C\ 1 

12:19:19.94 


/in /in 00 ""7 

49:40:22.7 


1.980 


23 


.00 


1 

Zt 





17 


n _i_ n on 
0.7o ± U.29 





.44 


_l_ 

zt 


n 
u 


1 1 

. 14 


HDr-MJVlD12 


12:37:19.80 


62:09:56.0 


2.648 


24. 


.on 


1 

Zt 


0. 


17 




1. 


02 


_l_ 
zt 


n 


zu 


tj tpi T7i D a \ r7 1 1 crt; 
HDr -BJVlZlloD 


TO '17 fi A O /I 

12: 67:04:. 64 


CO. 1 A A C O 

o2:14:4o.3 


2.211 


24 


.OZ 


1 

Zt 


0. 


20 







37 


_l_ 
zt 


n 
u. 


OQ 
.Uo 




1 0.Q7.0Q 1 C 

12:6 ( .26. 10 


o2:15:3o.U 


2.243 


23 


QP. 
. ;?o 


4- 
Zt 


0. 


1 A 
11 


n yic 1 n r\p 





.49 


-1- 

ZL 


n 


08 


J±L).b -BXloO 


T O O T On AT 

12:37:20.07 


/^O 1 O OO T 

02:12:22.7 


2.461 


24 


09 
OZ 


Zt 





17 


n no 1 n 1 r 

(J. 92 ± U.15 





74 


_l_ 

ZL 


n 
u. 


"1 9 
. 1Z 


Westpnal-JVLM47 


T /I TTCT r 7 0n 

14:17:57.39 


COOT C\ A cr 

52:31:04.5 


3.027 


24. 


.30 


± 


0. 


13 




1 


26 


± 


0. 


.10 


(^1422-073 


-\ A O 1 /I Z" 1 /II 

14:24:40.41 


00 cr cr /i cr cr 

22:55:45.5 


3.382 


24. 


.88 


=b 





.14 




1. 


05 


zb 


0. 


.17 




14-94-49 o°. 


99-o4-4n' ft 


2.229 


23 


.69 


zt 


Q 


07 




Q 


56 


zb 


0. 


.08 


Q1623-BX454 


16:25:51.42 


26:43:46.3 


2.422 


23. 


.89 


zt 


0. 


13 


0.50 ± 0.14 





22 


zb 


0. 


.06 


Q1623-BX663 


16:26:04.58 


26:48:00.2 


2.435 


24 


.14 


zt 


0. 


15 


1.02 ± 0.17 


0. 


24 


zb 


0. 


.07 


Q1623-BX747 


16:26:13.46 


26:45:53.2 


2.441 


22 


.55 


zt 





14 


0.42 ± 0.04 





12 


zb 


0. 


.03 


Q1700-MD157 


17:00:52.19 


64:15:29.3 


2.295 


24 


.35 


± 


0. 


13 


1.65 ± 0.28 





35 


zb 


0. 


.07 


Q1700-MD174 


17:00:54.54 


64:16:24.8 


2.347 


24 


.56 




0. 


17 


1.50 ± 0.28 





32 


zb 


0. 


,08 


Q2233-D3 


22:36:16.12 


13:55:19.2 


2.795 


24 


.08 







18 




1. 


08 


zb 


0. 


.16 


Q2233-MD21 


22:36:35.83 


13:55:42.0 


2.549 


25. 


.72 


± 


0. 


32 







85 


zb 


0. 


.21 


DSF2237A-D11 


22:40:02.99 


11:52:13.9 


2.959 


25 


.09 


zt 


0. 


26 







41 


zb 


0. 


.13 


DSF2237B-MD53 


22:39:28.67 


11:52:09.5 


2.292 


24 


.05 


± 


0. 


13 







.35 


zb 


0. 


.10 


Q2343-BX333 


23:46:21.51 


12:47:03.2 


2.397 


24 


.12 


± 


0. 


14 


0.65 ± 0.13 





30 


zb 


0. 


.08 


Q2346-BX445 


23:48:13.20 


00:25:15.8 


2.330 


23 


.81 


zb 





13 


1.00 ± 0.36 





39 


zb 


0. 


.09 



OBJECT TZ~J b K-K h [3.6 fan]" [4-5 A*m] a [5.8 /im] a [8.0 /<mf 



Q0000-C7 


2.17 


zb 


0. 


43 


4.50 


zb 


0.20 




































Q0000-C14 


4.43 


zb 


0. 


14 










































CDFB-D3 


4.97 


zb 


0. 


34 










































Q0100-BX172 










3.36 


zb 


0.22 


21.43 


zb 


0. 


10 


21 


.24 


zb 


0. 


10 


21 


.00 


zb 


0.23 


20.64 


zb 


0. 


10 


Q0142-BX195 










3.48 


zb 


0.18 


21.22 


zb 


0. 


10 


21 


.01 


zb 


0. 


10 


20 


.75 


zb 


0.24 


20.26 


zb 


0. 


10 


Q0142-BX256 










4.02 


zb 


0.13 


21.45 


zb 





10 


21 


.21 


zb 


0. 


18 


20 


.15 


zb 


0.34 










Q0201-OC12 


2.61 


zb 


0. 


40 


5.51 


zb 


0.20 




































Q0256-MD37 










3.68 


zb 


0.21 




































Q0933-MD38 










2.37 


zb 


0.20 




































Q1217-BX46 


2.36 


zb 


0. 


14 


3.56 


zb 


0.53 


21.65 


zb 


0. 


10 




























HDF-MMD12 


1.85 


zb 


0. 


27 


4.07 


zb 


0.21 


21.40 


zb 


0. 


07 


20 


.95 


zb 


0. 


07 


20 


.24 


zb 


0.07 


19.55 


zb 


0. 


07 


HDF-BMZ1156 


2.58 


zb 


0. 


20 


4.29 


zb 


0.17 


22.01 


zb 


0. 


07 


21 


.85 


zb 


0. 


07 


20 


.96 


zb 


0.07 


19.75 


zb 


0. 


07 


HDF-BMZ1384 


2.09 


zb 


0. 


16 


4.11 


zb 


0.14 


21.97 


zb 





07 


21 


.60 


zb 





07 


21 


.06 


zb 


0.07 


20.49 


zb 


0. 


07 


HDF-BX160 


2.09 


zb 


0. 


18 


3.97 


zb 


0.20 


21.39 


zb 


0. 


07 


21 


22 


zb 


0. 


07 


20 


.90 


zb 


0.07 


20.85 


zb 


0. 


07 


Wcstphal-MM47 
















21.83 


zb 





10 


21 


.84 


zb 


0. 


10 


















Q1422-C73 










3.82 


zb 


0.21 




































Q1422-MD109 










4.48 


zb 


0.10 




































Q1623-BX454 
















22.87 


zb 





15 


22 


.64 


zb 





15 


















Q1623-BX663 


1.63 


zb 


0. 


26 


4.22 


zb 


0.14 


20.66 


zb 


0. 


10 


19 


.86 


zb 


0. 


TO 


19 


.03 


zb 


0.10 


18.13 


zb 


0. 


10 


Q1623-BX747 


1.51 


zb 


0. 


11 


2.19 


zb 


0.23 


22.33 


zb 


0. 


16 




























Q1700-MD157 


2.24 


zb 


0. 


10 


4.12 


zb 


0.19 










21 


.61 


zb 


0. 


10 










21.25 


zb 


0. 


11 


Q1700-MD174 










4.66 


zb 


0.15 










20 


.41 


zb 


0. 


10 










18.26 


zb 


0. 


10 


Q2233-D3 










4.53 


zb 


0.19 




































Q2233-MD21 










4.37 


zb 


0.31 


22.89 


zb 


0. 


11 




























DSF2237A-D11 










3.00 


zb 


0.50 




































DSF2237B-MD53 










4.68 


zb 


0.14 




































Q2343-BX333 


1.93 


zb 


0. 


14 


4.21 


zb 


0.16 










21 


.01 


zb 


0. 


10 










20.47 


zb 


0. 


10 


Q2346-BX445 










4.34 


zb 


0.13 





































a U , G, 7Z. and IR.AC magnitudes arc on the AB system, 
k 71 — J and 1Z — K colors are AB — Vega. 



HAINLINE ET AL. 



TABLE 2 

Best Fit Parameters, AGNs with IRAC Coverage 



OBJECT 


log [Age] 
Myr 


log(M») 
M 


E(B-V) SF 


SFR 

M Q yr- 1 


E(B - 


- V)agn 




O.Z1 It U.lo 


10 


.77 ± 0. 


.09 


U.14U it U.UoU 


Q7 _|_ 


z.DUU 


± 


2. 


.433 


Q0142-BX195 


2.76 ± 0.31 


10 


.76 ± 0. 


15 


0.260 ± 0.053 


101 ± 62 


4.600 


± 


2. 


.066 


Q0142-BX256 


3.44 ± 0.08 


10 


.96 ± 


.06 


0.180 ± 0.021 


33 ± 6 


5.600 


± 


1. 


010 


HDF-MMD12 


2.31 ± 0.11 


10 


.60 ± 


12 


0.400 ± 0.033 


196 ± 49 


2.000 


± 





.614 


HDF-BMZ1156 


3.34 ± 0.10 


10 


.64 ± 0. 


.05 


0.240 ± 0.024 


20 ± 3 


6.200 


± 


0. 


646 


HDF-BMZ1384 


3.44 ± 0.00 


10 


.69 ± 


.07 


0.140 ± 0.018 


18 ± 3 


2.800 


± 


1. 


.254 


HDF-BX160 


3.11 ± 0.21 


10 


.91 ± 0. 


.08 


0.260 ± 0.029 


64 ± 15 


6.400 


± 


1. 


.231 


Q1623-BX663 


2.76 ± 0.81 


10 


.26 ± 


.85 


0.180 ± 0.020 


32 ± 4 


1.000 


± 


0. 


.078 


Q1700-MD157 


3.16 ± 0.17 


10 


.74 ± 0. 


.09 


0.260 ± 0.025 


38 ± 8 


7.400 


± 


0. 


.616 


Q1700-MD174 


3.44 ± 0.00 


11 


.15 ± 0. 


17 


0.300 ± 0.042 


51 ± 16 


3.600 


± 


1 


.361 


Q2343-BX333 


3.41 ± 0.00 


10 


.96 ± 


09 


0.200 ± 0.019 


35 ± 7 


1.000 


± 


2. 


.461 




io gV iv AGiV } 


SAGN, l M m b 


lOg^bol.AGN ) 


Y 2 v 2 c 










Q0100-BX172 


-1.00 ± 0.55 




0.02 




10.82 ± 0.57 


0.26, 0.09 










Q0142-BX195 


-1.00 ± 0.27 




<0.01 




11.26 ± 0.35 


1.03^ 0.34 










Q0142-BX256 


0.18 ± 0.00 




0.01 




11.95 ± 0.09 


6.30, 3.15 










HDF-MMD12 


-1.00 ± 0.00 




0.16 




11.55 ± 0.08 


4.63, 1.54 










HDF-BMZ1156 


0.15 ± 0.07 




0.01 




11.70 ± 0.08 


19.06, 6.35 










HDF-BMZ1384 


-0.52 ± 0.11 




0.04 




10.98 ± 0.15 


45.06, 11.26 










HDF-BX160 


-1.00 ± 0.00 




<0.01 




11.06 ± 0.11 


15.39, 3.85 










Q1623-BX663 


0.18 ± 0.02 




0.79 




11.94 ± 0.04 


4.91, 1.23 










Q1700-MD157 


-1.00 ± 0.20 




<0.01 




10.84 ± 0.18 


22.97, 11.48 










Q1700-MD174 


0.15 ± 0.02 




0.09 




12.11 ± 0.14 


13.33, 13.33 










Q2343-BX333 


-1.00 ± 0.00 




0.11 




10.80 ± 0.10 


14.11, 7.06 











a N AGN = L 6oi /(SFR x 10 10 ) 



The fraction of AGN emission at rcst-framc 1 fim. 
c Raw and reduced \ 2 values for each best-fit model. 



TABLE 3 

Best Fit Parameters, AGNs with IRAC Coverage, SPS-only, through K-band 



OBJECT 


log [Age] 


log(M») 


E(B-V)sf SFR 


Y 2 Y 2 , a 




Myr 


M© 


M yr- 1 





Q0100-BX172 


3.21 


± 


0. 


.23 


10 


.77 


± 


0. 


20 





.140 


± 





.032 


37 


± 


12 


0.19, 


0.19 


Q0142-BX195 


2.66 


± 





32 


10 


.68 


± 


0. 


17 





.265 


± 


0. 


056 


106 


± 


50 


0.77, 


0.77 


Q0142-BX256 


3.44 


± 


0. 


00 


11 


.09 


± 


0. 


11 


0. 


.210 


± 





019 


45 


± 


9 


1.81, 


1.81 


HDF-MMD12 


1.86 


± 


0. 


.72 


10 


.56 


± 


0. 


27 


0. 


.485 


± 


0. 


.112 


505 


± 


757 


0.16, 


0.16 


HDF-BMZ1156 


3.44 


± 





00 


10 


.95 


± 


0. 


26 





.265 


± 





.024 


32 


± 


10 


8.44, 


8.44 


HDF-BMZ1384 


3.44 


± 


0. 


00 


10 


.91 


± 


0. 


14 


0. 


.175 


± 


0. 


.016 


29 


± 


6 


34.12, 


17.06 


HDF-BX160 


3.34 


± 





13 


11 


.17 


± 


0. 


26 


0. 


.255 


± 





028 


68 


± 


24 


12.73, 


6.37 


Q1623-BX663 


3.41 


± 





00 


11 


.06 


± 


0. 


66 





.220 


± 


0. 


.019 


44 


± 


11 


9.83, 


4.91 


Q1700-MD157 


3.44 


± 


0. 


00 


11 


.10 


± 


0. 


25 


0. 


.265 


± 





.021 


45 


± 


12 


16.13, 


8.07 


Q1700-MD174 


3.44 


± 





00 


11 


.14 


± 


0. 


12 


0. 


.295 


± 


0. 


.022 


51 


± 


14 


13.31, 


13.31 


Q2343-BX333 


3.41 


± 


0. 


.00 


11 


.02 


± 


0. 


11 


0. 


.210 


± 





.018 


40 


± 


9 


14.91, 


7.45 



a Raw and reduced \ 2 values for each best-fit model. 



UV-SELECTED AGN STELLAR MODELING 



TABLE 4 

AGN Correction Factors 



OBJECT 


Mass a 


SFR a 


Agc a 


E(B 


— 


V) S F a 


Q0100 


- BX172 


i n -t- n ^ 
±.u zn u.o 


1.0 


± 0.4 


1.0 


± 0.6 


1.0 


± 


0.3 


Q0142 


- BX195 


1.2 ± 0.6 


1.1 


± 0.7 


1.3 


± 1.0 


1.0 


± 


0.3 


Q0142 


- BX256 


0.7 ± 0.2 


0.8 


± 0.3 


0.8 


± 0.2 


0.9 


± 


0.1 


HDF - 


MMD12 


0.9 ± 0.7 


0.2 


± 1.1 


2.1 


± 1.4 


0.8 


± 


0.2 


HDF - 


BMZ1156 


0.5 ± 0.3 


0.6 


± 0.3 


0.8 


± 0.2 


0.9 


± 


0.1 


HDF - 


BMZ1384 


0.6 ± 0.3 


0.6 


± 0.4 


1.0 


± 0.0 


0.8 


± 


0.2 


HDF - 


BX160 


0.6 ± 0.4 


0.8 


± 0.4 


0.8 


± 0.5 


1.0 


± 


0.2 


Q1623 


- BX663 


0.3 ± 1.0 


0.8 


± 0.3 


0.4 


± 0.9 


0.8 


± 


0.1 


Q1700 


-MD157 


0.5 ± 0.3 


0.8 


± 0.3 


0.6 


± 0.4 


1.0 


± 


0.1 


Q1700 


- MD174 


0.9 ± 0.5 


0.9 


± 0.4 


1.0 


± 0.0 


1.0 


± 


0.2 


Q2343 


- BX333 


1.0 ± 0.3 


1.0 


± 0.3 


1.0 


± 0.0 


1.0 


± 


0.1 


average 13 


0.7 ± 0.3 


0.8 


± 0.2 


1.0 ± 0.4 


0.92 ± 0.07 



a Correction values represent the derived parameters from the AGN+SF 
modeling fit divided by the SF modeling fit through the K band data, 
k Uncertainties for the average correction factors represent the standard 
deviation of the distributions. 



TABLE 5 



Best Fit Parameters, AGNs without IRAC Coverage 




log [Age] 


log [Age 


] corr 


log(M„) 


log(M„) corr 








Myr 


Myr 










O0000-C7 


3.26 


± 0.00 


3.25 ± 


0.19 


10 35 + 14 


10.21 ± 0.29 






Q0000-C14 


3.30 


± 0.00 


3.30 ± 


0.19 


11.43 ± 0.12 


11.30 ± 0.27 






CDFB-D3 


3.36 


± 0.00 


3.36 ± 


0.19 


11.15 ± 0.21 


11.02 ± 0.36 






Q0201-OC12 


3.41 


± 0.00 


3.41 ± 


0.19 


11.53 ± 0.19 


11.39 ± 0.34 






Q0256-MD37 


3.06 


± 0.29 


3.05 ± 


0.35 


10.78 ± 0.21 


10.65 ± 0.36 






Q0933-MD38 


2.96 


± 0.26 


2.95 ± 


0.33 


10.52 ± 0.25 


10.39 ± 0.41 






Q1217-BX46 


2.86 


± 0.14 


2.85 ± 


0.24 


10.49 ± 0.08 


10.35 ± 0.25 






WEST-MM47 


2.06 


± 0.35 


2.05 ± 


0.41 


10.41 ± 0.13 


10.27 ± 0.28 






Q1422-C73 


3.26 


± 0.00 


3.25 ± 


0.19 


10.82 ± 0.15 


10.69 ± 0.30 






Q1422-MD109 


3.01 


± 0.27 


3.00 ± 


0.34 


11.29 ± 0.14 


11.16 ± 0.29 






Q1623-BX454 


2.76 


± 0.31 


2.75 ± 


0.36 


10.03 ± 0.19 


9.90 ± 0.33 






Q1623-BX747 


2.81 


± 0.12 


2.80 ± 


0.23 


10.32 ± 0.09 


10.19 ± 0.25 






Q2233-D3 


2.71 


± 0.63 


2.70 ± 


0.66 


11.36 ± 0.25 


11.22 ± 0.41 






Q2233-MD21 


2.66 


± 0.71 


2.65 ± 


0.74 


10.20 ± 0.19 


10.07 ± 0.34 






DSF2237A-D11 


3.32 


± 0.16 


3.32 ± 


0.25 


10.20 ± 0.20 


10.06 ± 0.35 






DSF2237B-MD53 


3.44 


± 0.00 


3.44 ± 


0.19 


11.37 ± 0.10 


11.24 ± 0.26 






Q2346-BX445 


3.44 


± 0.00 


3.44 ± 


0.19 


11.33 ± 0.10 


11.19 ± 0.26 






OBJECT 


E(B 


-V) SF 


E(B - V) 


SF,corr 


SFR 


SFRcorr 


2 2 a,b 




M 


M 


M yr- 1 


M Q yr" 1 






Q0000-C7 


0.015 


± 0.017 


0.014 ± 


1.133 


12 ± 4 


10 ± 5 


151.65, 


151.65 


Q0000-C14 


0.315 


± 0.024 


0.289 ± 


0.076 


135 ± 37 


104 ± 55 


2.56, 




CDFB-D3 


0.295 


± 0.044 


0.271 ± 


0.149 


62 ± 30 


48 ± 36 


8.20, 




Q0201-OC12 


0.460 


± 0.032 


0.422 ± 


0.070 


129 ± 57 


99 ± 68 


7.44, 


7.44 


Q0256-MD37 


0.225 


± 0.052 


0.206 ± 


0.231 


53 ± 31 


41 ± 35 


0.00, 




Q0933-MD38 


0.000 


± 0.000 


0.000 ± 


0.000 


37 ± 2 


28 ± 11 


0.94, 




Q1217-BX46 


0.255 


± 0.026 


0.234 ± 


0.102 


43 ± 10 


33 ± 16 


3.88, 


1.94 


WEST-MM47 


0.355 


± 0.026 


0.326 ± 


0.073 


225 ± 34 


174 ± 75 


6.84, 


6.84 


Q1422-C73 


0.185 


± 0.031 


0.170 ± 


0.168 


37 ± 13 


28 ± 17 


2.84, 




Q1422-MD109 


0.380 


± 0.037 


0.348 ± 


0.097 


194 ± 70 


150 ± 91 


10.67, 


10.67 


Q1623-BX454 


0.090 


± 0.045 


0.083 ± 


0.500 


19 ± 5 


15 ± 7 


3.68, 


1.84 


Q1623-BX747 


0.030 


± 0.010 


0.028 ± 


0.333 


33 ± 5 


25 ± 11 


43.23, 


14.41 


Q2233-D3 


0.440 


± 0.096 


0.403 ± 


0.218 


446 ± 303 


344 ±331 


0.00, 




Q2233-MD21 


0.380 


± 0.046 


0.348 ± 


0.121 


35 ± 7 


27 ± 13 


1.47, 


1.47 


DSF2237A-D11 


0.065 


± 0.050 


0.060 ± 


0.769 


8 ± 6 


6 ± 6 


0.00, 




DSF2237B-MD53 


0.310 


± 0.022 


0.284 ± 


0.071 


85 ± 21 


66 ± 33 


8.21, 




Q2346-BX445 


0.265 


± 0.021 


0.243 ± 


0.079 


77 ± 19 


60 ± 30 


2.26, 


2.26 



Raw and reduced x 2 values for each best-fit model. 

Objects with only three photometric data points have reduced \ 2 values given by a dash. 
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TABLE 6 

UV Emission Features for Stellar Mass Subsamples 





log(M»/Af ) > 10.7 


log(M*/M©) < 10.7 


w LyQ a 

WNV,1240 a 
^CIV,1549 a 
M / HoII,1640 a 


56±15 
4.6±1.4 
17.6±3.3 
10.2±2.4 


58±18 
3.7±1.6 
7.4±3.3 
5.0±1.1 



a Rcst-framc EW in A, measured from the composite spectra. 
Uncertainties are calculated as described in j|5.4| 



